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ABSTRACT 

Triggered by the earlier surprising detection of pulsed hard X-ray emission from IE 1841-045, we investigated 
the time-averaged high-energy spectral characteristics of the established Anomalous X-ray Pulsars, 1RXS J 1708 - 
4009, 4U 0142 + 61, IE 2259 + 586 and IE 1048.1-5937, all with persistent X-ray emission. We report on 
the discovery of hard spectral tails for energies above 10 keV in the total and pulsed spectra of AXPs 1RXS 
J1708 -4009, 4U 0142 + 61 and IE 2259 + 586 using RXTE PCA (2-60 keV) / HEXTE (15-250 keV) data and 
INTEGRAL IBIS ISGRI (20-300 keV) data. IE 1048.1-5937 appeared to be too weak to be detected with the 
presently available exposure. Furthermore, improved spectral information for IE 1841-045 is presented. The 
pulsed and total spectra measured above 10 keV have power-law shapes and there is so far no significant evidence 
for spectral breaks or bends up to ~ 150 keV. The pulsed spectra above 10 keV are exceptionally hard with indices 
measured for four AXPs approximately in the range -1.0 - 1.0. The indices measured for these pulsed spectra 
below 10 keV were in the range 2.0 - 4.3, indicating the very drastic spectral changes in a narrow energy interval 
around 10 keV, where clear minima are found in luminosity. The best fit power-law models to the total spectra 
between ~ 10 and 150 keV are significantly softer, with indices measured for three AXPs, IE 1841-045, 1RXS 
J 1708 -4009 and 4U 0142 + 61 in the range 1.0- 1.4. For the latter AXPs the pulsed fractions are consistent with 
100% around 1 00 keV, but are different at 10keV:4U 0142 + 61 - 10%, IE 1841-045 ~25%, and 1RXS J1708- 
4009 consistent with 100%. The luminosities of these total and pulsed spectral tails (10- 150 keV) largely exceed 
the total available spin-down powers by factors ranging from ~ 100 to ~ 600. This shows that also these new 
hard-X-ray components cannot be powered by rotational energy loss. We reanalyzed archival CGRO COMPTEL 
(0.75-30 MeV) data to search for signatures from our set of AXPs. No detections can be claimed, but the obtained 
upper-limits to soft gamma-ray emission in this MeV range indicates for 1RXS J1708 -4009, 4U 0142 + 61 and 
IE 1841-045 that strong breaks or bends must occur somewhere between ^150 keV and 750 keV. We discuss 
predictions from first attempts to model our hard X-ray / soft gamma-ray spectra in the context of the magnetar 
model. Our spectral results cannot yet discriminate between the different proposed scenarios. Particularly, more 
constraining information is required on the detailed shape of the spectra between ~ 150 keV (our highest data 
points from INTEGRAL) and 750 keV, from where we report the COMPTEL upper limits. 

Subject headings: pulsars: individual (1RXS J1708-4009; 4U 0142 + 61; IE 2259 + 586;lE 1048.1 -5937;1E 
1841-045), X-rays: stars 



1. INTRODUCTION 

Anomalous X-ray pulsars (AXP) belong to a class of rare 
objects closely concentrated along the Galactic plane, known 
to emit pulsed X-rays for energies below 10 keV with pulse pe- 
riods in the range ~ 6 - 12 s and characteristic spin-down time 
scales of~10 3 -10 5 year. Two, probably three, members of this 
class are embedded in a shell-like supernova remnant (SNR; 
see e.g. Gregory & Fahlman 1980; Kriss et al. 1985). The fact 
that the observed X-ray luminosity is much larger than the spin- 
down power excludes an interpretation in which the (pulsed) X- 
ray emission originates from a spin-down powered pulsar. On 
the other hand, the steady spin-down without a doppler modu- 
lated signature and the lack of bright optical counterparts make 
an X-ray binary interpretation, in which mass transfer (accre- 
tion) powers the high-energy emission of the system, very un- 
likely (see Mereghetti et al. 2002, for a review on AXPs). Cur- 
rently, models based on the decay of very strong magnetic fields 
(10 14 -10 15 GauB) - so called "magnetar" models (Thompson & 
Duncan 1996) - seem to explain the observed high-energy char- 
acteristics of AXPs at a satisfactory level. For instance, the re- 
cently detected bursts from the AXPs IE 1048.1 - 5937 (Gavriil 
et al. 2002b) and IE 2259 + 586 (Kaspi et al. 2003) mimic the 
bursting behaviour of Soft Gamma-ray Repeaters (SGR) for 
which the magnetar model was initially developed. Also, the 
"glitch" phenomenon detected in the spin-down of some AXP 
members fits in this model (Kaspi et al. 2000, 2003; Morii et 



al. 2005). These properties provide strong evidence that both 
AXPs and SGRs are members of the same source class. 

The X-ray spectra of AXPs in the 0.5-10 keV band are very 
soft and can best be described by a black body plus a power-law 
model. The softness of the spectra below 10 keV (power-law 
indices 2 < T < 4, with F 1 oc E~ T ) predicts non-detections for 
energies above 10 keV and thus explains the initial ignorance 
for studying the spectral properties of AXPs at energies above 
lOkeV. 

It was a great surprise that the high-resolution imaging in- 
strument IBIS/ISGRI aboard ESAs INTEGRAL satellite mea- 
sured hard X-rays from the direction of three AXPs. Firstly, 
Molkov et al. (2004) reported the discovery of an INTEGRAL 
source at the position of AXP IE 1841-045 in SNR Kes 73 
for energies up to 120 keV (60-120 keV: 7.5 a). This was fol- 
lowed up by Kuiper et al. (2004), who analysed archival RXTE 
PCA and HEXTE data from monitoring observations spread 
over four years, to prove that the hard X-ray emission comes 
from the AXP and not from the SNR. They discovered non- 
thermal pulsed hard X-ray / soft 7-ray emission up to ~ 150 
keV with a spectrum with power-law photon index of ~ 0.94. 

Secondly, Revnivtsev et al. (2004) published the INTEGRAL 
detection ( 1 8-60 keV: 6.5 cr) of AXP 1 RXS J 1 70849.0 - 4009 1 
(we use throughout this paper: 1RXS J1708-4009), also in 
a spatial analysis of ISGRI data. Similarly, den Hartog et al. 
(2004) reported the detection of AXP 4U 0142 + 61 (50-100 



1 



2 



Kuiper et al. 



keV: 6.1 a) which exhibited a very hard total spectrum above 20 
keV. These three AXPs reach about the same flux level around 
100 keV. 

In this paper we present the results from follow-up studies 
using archival RXTE PC A and HEXTE data of 1RXS J 1708- 
4009, 4U 0142 + 61, IE 2259 + 586 and IE 1048.1 -5937 aimed 
at studying their timing and spectral characteristics above 10 
keV. We also revisit the pulsed high-energy emission properties 
of IE 1841-045 above 10 keV using more RXTE PCA data and 
applying now more-optimized event selection criteria. Further- 
more, we explored the INTEGRAL database using both public, 
private and core program data to derive the total X-ray emis- 
sion spectra of the five AXPs mentioned above, and reanalyze 
archival data from COMPTEL (Schonfelder et al. 1993), de- 
livering constraining upper limits in the soft gamma-ray band 
0.75-30 MeV. Finally, initial results will be shown from IBIS 
ISGRI timing analysis studies of 1RXS J1708-4009 and IE 
1841-045. 

2. INSTRUMENTS AND OBSERVATIONS 
2. 1 . Rossi X-ray Timing Explorer 

In this study extensive use is made of data from monitoring 
observations of AXPs with the two non-imaging X-ray instru- 
ments aboard RXTE, the Proportional Counter Array (PCA; 2- 
60 ke V) and the High Energy X-ray Timing Experiment (HEXTE; 
15-250 keV). The PCA (Jahodaetal. 1996) consists offivecol- 
limated xenon proportional counter units (PCUs) with a total 
effective area of - 6500 cm 2 over a - 1° (FWHM) field of 
view. Each PCU has a front Propane anti-coincidence layer and 
three Xenon layers which provide the basic scientific data, and 
is sensitive to photons with energies in the range 2-60 keV. The 
energy resolution is about 18% at 6 keV. 

The HEXTE instrument (Rothschild et al. 1998) consists of 
two independent detector clusters, each containing four Na(Tl)/ 
CsI(Na) scintillation detectors. The HEXTE detectors are me- 
chanically collimated to a ~ 1° (FWHM) field of view and 
cover the 15-250 keV energy range with an energy resolution 
of ~ 15% at 60 keV. The collecting area is 1400 cm 2 taking 
into account the loss of the spectral capabilities of one of the 
detectors. The maximum time resolution of the tagged events 
is 7.6/is. In its default operation mode the field of view of each 
cluster is switched on and off source to provide instantaneous 
background measurements. Due to the co-alignment of HEXTE 
and the PCA, they simultaneously observe the sources. Table 
1 lists the publicly available RXTE observations used in this 
study. In the fourth column the PCU unit-2 screened (see Sect 
3.1) exposure is given. A typical observation consists of sev- 
eral sub-observations spaced more or less uniformly between 
the start and end date of the observation. 

2.2. INTEGRAL 

The INTEGRAL spacecraft (Winkler et al. 2003), launched 
17 October 2002, carries two main 7-ray instruments: a high- 
angular-resolution imager IBIS (Ubertini et al. 2003) and a high- 
energy-resolution spectrometer SPI (Vedrenne et al. 2003). These 
instruments make use of coded aperture masks enabling image 
reconstruction in the hard X-ray/soft 7-ray band. 

In our study, guided by sensitivity considerations, we only 
used data recorded by the INTEGRAL Soft Gamma-Ray Im- 
ager ISGRI (Lebrun et al. 2003), the upper detector system of 
IBIS, sensitive to photons with energies in the range ~20 keV 
- 1 MeV. With an angular resolution of about 12' and a source 



Table 1 

List of RXTE observations of 1RXS J1708 - 4009, 4U 
0142 + 61, IE 2259 + 586 and IE 1048. 1 - 5937 used in this 
study. For observations used in the reanalysis of IE 
1841-045 see Sect. 8.1 and Kuiper et al. (2004). 



Obs. 


Begin/End Date 


17 t 

Exp. 1 


id. 


(dd/mm/yyyy) 


(ks) 


1RXS Jl 70S -4009 






30125 


12-01-1998 


08-01-1999 


59.896 


40083 


06-02-1999 


1 1 fiQ onnn 
1 l-Uo-zUUU 


52.568 


50082 


21-04-2000 


12-UJ-2UU1 


34.576 


60412 


20-05-2001 


T2 KVY1 

2J-UJ-ZUU1 


9.928 


60069 


06-05-2001 


Zxj-UZ-ZOVZ 


24.544 


70094 


02-04-2002 


Zu-Uj-ZUUj 


55 600 


80098 


16-04-2003 


2o- 1 U-2UU J 


73.368 


All 


12-01-1998 


Zo- 1 U-ZOUS 


310 480 


41] 0142 + 61 






10193 


28-03-1996 


29-03-1996 


37.728 


10185 


28-03-1996 


28-03-1996 


16.288 


20146 


24-11-1996 


13-12-1997 


10.600 


30110 


21-03-1998 


21-03-1998 


15.408 


50082 


07-03-2000 


1 n no onn 1 
1U-U2-ZU01 


34.240 


60069 


18-03-2001 


no ni irwi 
08-01-2002 


42.272 


70094 


06-03-2002 


O^ 1 O 

Zo-12-ZUUZ 


82.464 


80098 


28-03-2003 


1 q no 

1 ©-U9-ZUU J 


38.664 


80099 


03-09-2003 


no no iivn 


29.216 


All 


28-03-1996 


1 q no onna 
1 8-09-2003 


306 880 

JV/VJ. OOV7 


IE 2259 + 586 






10192 


29-09-1996 


30-09-1996 


75.936 


20145 


25-02-1997 


25-03-1997 


101.024 


20146 


24-11-1996 


13-12-1997 


9.784 


30126 


13-08-1998 


02-12-1998 


103.136 


40083 


17-01-1999 


01-03-2001 


A Q -I A A 

45.344 


40082 


26-01-2000 


27-03-2000 


87.120 


50082 


10-03-2000 


02-03-2001 


50.280 


60069 


16-04-2001 


07-02-2002 


49.592 


70094 


22-03-2002 


15-02-2003 


149.840 


80098 


15-03-2003 


28-10-2003 


71.960 


All 


29-09-1996 


28-10-2003 


747.016 


IE 1048.1-5937 






10192 


29-07-1996 


30-07-1996 


72.344 


20146 


24-11-1996 


13-12-1997 


12.616 


40083 


23-01-1999 


10-02-2000 


49.184 


50082 


11-03-2000 


09-02-2001 


40.992 


60069 


06-03-2001 


25-02-2002 


139.776 


70094 


12-03-2002 


25-02-2003 


133.224 


80098 


12-03-2003 


24-02-2004 


136.376 


All 


29-07-1996 


24-02-2004 


584.512 


tPCU-2 exposure after screening 



location accuracy of better than 1' (for a > 10cr source) this 
instrument is able to locate and separate high-energy sources 
in crowded fields within its 19° x 19° field of view (50% par- 
tially coded) with an unprecedented sensitivity (~ 960 cm 2 at 
50 keV). Its energy resolution of about 7% at 100 keV is am- 
ply sufficient to determine the (continuum) spectral properties 
of hard X-ray sources in the ~ 20 - 300 keV energy band. 

The timing accuracy of the ISGRI time stamps recorded on 
board is about 61^is. The time alignment between INTEGRAL 
and RXTE is better than ~ 25^s, verified using data from si- 
multaneous RXTE and INTEGRAL observations of the accretion- 
powered millisecond pulsar IGR J0029 1+5934 (Falanga et al. 
2005); for a calibration on the Crab pulsar, see also Kuiper et 
al. (2003). Given the fact that the accuracy of the RXTE clock 
in absolute time is about 2^s (Rots et al. 2004), this implies that 
the INTEGRAL absolute timing is better than ~ 27 fis. Data 
from regular INTEGRAL Crab monitoring observations show 
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Table 2 

List of INTEGRAL observations, sorted on INTEGRAL 

ORBITAL REVOLUTIONS (REV.), OF THE AXPS STUDIED IN THIS 

work. For more details on the executed INTEGRAL 

OBSERVATIONS, SEE HTTP://INTEGRAL.ESAC.ESA.INT/ 



Rev. Rev. Begin/End Date 

begin end (dd/mm/yyyy) 

Imaging analysis 



1RXS J 1708 -4009 

36 106 29-01-2003 29-08-2003 



41] 0142 + 61 






47 


92 


03-03-2003 


16-07-2003 


142 


148 


12-12-2003 


01-01-2004 


153 


153 


14-01-2004 


15-01-2004 


161 


162 


07-02-2004 


12-02-2004 


177 


234 


26-03-2004 


12-09-2004 


238 


238 


25-09-2004 


27-09-2004 


261 


266 


02-12-2004 


19-12-2004 


268 


269 


24-12-2004 


28-12-2004 


IE 1841-045 






49 


253 


10-03-2003 


08-11-2004 


IE 2259 + 586 






142 


148 


12-12-2003 


01-01-2004 


161 


162 


07-02-2004 


12-02-2004 



IE 1048.1-5937 

36 217 29-01-2003 24-07-2004 



Timing analysis 

1RXS Jl 708 -4009 = 
36 120 29-01-2003 10-10-2003 

IE 1841-045 

49 123 10-03-2003 18-10-2003 



that the clock behaviour is stable, making deep timing studies 
of weak pulsars possible. 

In its default operation mode INTEGRAL observes the sky in 
a dither pattern with 2° steps, which could be rectangular e.g. 
a 5 x 5 dither pattern with 25 grid points, or hexagonal with 
7 grid points (target in the middle). Typical integration times 
for each grid point (pointing/sub-observation) are in the range 
1800 - 3600 seconds. This strategy drives the structure of the 
INTEGRAL data archive which is organised in so-called sci- 
ence windows (Sew) per INTEGRAL orbital revolution (last- 
ing for about 3 days) containing the data from all instruments 
for a given pointing. Most of the INTEGRAL data reduction 
in this study was performed with the Offline Scientific Analysis 
(OS A) version 4. 1 distributed by the INTEGRAL Science Data 
Centre (ISDC; see e.g. Courvoisier et al. 2003). 

Table 2 lists the INTEGRAL orbital revolution identifiers 
with corresponding start/end dates of the observations used in 
the imaging/spectral analyses and timing analyses of the se- 
lected sample of persistent AXPs. 

2.3. COMPTEL 

The Compton telescope COMPTEL aboard the Compton Gan 
Ray Observatory (CGRO, 1991 - 2000) was sensitive to 7-ray 
photons between 0.75 and 30 MeV, thereby covering the harder 
7-ray band adjacent to the INTEGRAL one. The very hard 
spectra that we measured with IBIS ISGRI for some AXPs war- 
ranted us to revisit the COMPTEL data archive to search for 
signals from AXPs. COMPTEL has an energy-dependent en- 
ergy and angular resolution of 5% - 8% (FWHM) and l.°7 - 
4.° 4 (FWHM), respectively, and a wide circular field of view 
covering ~1 steradian. Imaging in its large field of view is pos- 



sible with a location accuracy (flux dependent) of the order of 
0.°5 - 2°. For details on the experiment see Schonfelder et al. 
(1993). 

3. ANALYSIS METHODS 

3.1. RXTE PCA/HEXTE timing 

The PCA data from the observations listed in Table 1 have all 
been collected in Goodxenon or GoodxenonwithPropane mode 
allowing high-time-resolution (0.9/xs) studies in 256 spectral 
channels. Because we are mainly interested in the medium/hard 
(> 2 keV) X-ray timing properties of the selected sample of 
AXPs, we ignored the events triggered in the Propane layers 
of each PCU. Furthermore, contrary to the work presented in 
Kuiper et al. (2004), we now used data from all three xenon 
layers of each PCU, namely, employing data from the (deeper) 
middle and lower xenon layers considerably improves the signal- 
to-noise ratio for energies above ^10 keV. This allows us to 
better characterize the hard (> 10 keV) X-ray properties of 
AXPs. 

The number of active PCUs at any time was changing, there- 
fore, we treated the five PCUs constituting the PCA separately. 
Good time intervals have been determined for each PCU by 
including only time periods when the PCU in question is on, 
and during which the pointing direction is within 0.°05 from 
the target, the elevation angle above Earth's horizon is greater 
than 5°, a time delay of 30 minutes since the peak of a South- 
Atlantic-Anomaly passage holds, and a low background level 
due to contaminating electrons is observed. These good time 
intervals have subsequently been applied in the screening pro- 
cess to the data streams from each of the PCUs (e.g. see Table 
1 for the resulting screened exposure of PCU-2 per observation 
run). 

Next, for each sub-observation the arrival times of the se- 
lected events (for each PCU unit) have been converted to arrival 
times at the solar system barycenter (in TDB (=barycentric dy- 
namical time) time scale; DE200 solar system ephemeris) using 
the instantaneous spacecraft position and the celestial positions 
of the selected sample of AXPs. In this work we used for the 
AXP positions those obtained by the Chandra X-ray observa- 
tory with a typical position accuracy of about 0"5 (for 1RXS 
J1708-4009, see Israel et al. (2003); for4U 0142 + 61, see Pa- 
tel et al. (2003); for IE 2259 + 586, see Patel et al. (2001); for 
IE 1048. 1 - 5937, see Wang & Chakrabarty (2002); Israel et al. 
(2002) ; and finally for IE 1841-045, see Wachter et al. (2004)). 

These bary centered arrival times have been folded with avail- 
able phase connected timing solutions (see for details the rele- 
vant section on the timing characteristics for each AXP) us- 
ing only the first three frequency coefficients (frequency, first 
and second frequency time derivatives at a certain epoch) to 
obtain pulse phase distributions for selected energy windows. 
Combining now the phase distributions from the various PCUs 
for energies between ~ 2 and 10 keV we obtained for each 
- sub-observation pulse profiles, which deviate from uniformity 
at significances well above 5ct for each of the AXPs in our 
sample. For the calculation of these significances we applied 
in this work the bin-free Z\ statistic (Buccheri et al. 1983), 
which behaves as a \ 2 distribution for 2n degrees of freedom 
(n=number of harmonics). However, phase shifts between the 
sub-observations made a direct combination of the pulse pro- 
files to obtain (very) high-statistics time-averaged pulse profiles 
impossible. Therefore we correlated the pulse phase distribu- 
tion of each sub-observation with a chosen initial template and 
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applied the measured phase shifts to obtain an aligned com- 
bination with much higher statistics. The correlation is then 
once repeated with, instead of the initially chosen template, the 
aligned combination from the first correlation to obtain the fi- 
nal summed profile (see e.g. de Plaa et al. 2003, for a similar 
iterative method applied for PSR B0540-69). 

The net result is a set of aligned high-statistics pulse profiles 
in 256 energy channels for each of the AXPs in our sample. It 
should be noted that these profiles are time-averaged profiles 
ignoring possible temporal variations in shape and/or pulsed 
fraction (see e.g. Rea et al. (2005) for 1RXS J1708 - 4009; Rea 
et al. (2006) for 4U 0142 + 61; and Gavriil & Kaspi (2004); 
Tiengo et al. (2005) for IE 1048.1 - 5937). 

HEXTE operated in its default rocking mode during the ob- 
servations listed in Table 1, allowing the collection of real-time 
background data from two independent positions ±1.°5 to ei- 
ther side of the on-source position. For the timing analysis 
we selected only the on-source data. Good time intervals have 
been determined using similar screening filters as used in the 
case of the PCA. The selected on-source HEXTE event times 
have subsequently been barycentered and folded according to 
the available ephemerides (see PCA part at the beginning of 
this section) using again only the first three frequency coeffi- 
cients. Applying, for each AXP in our sample, the phase shifts 
as derived from the contemporaneous PCA measurements to 
the HEXTE phase distributions of each sub-observation we ob- 
tained the time averaged HEXTE pulse phase distributions in 
256 spectral channels (15 - 250 keV) for the combination of 
observations listed in Table 1 . 

3.2. RXTE PCA/HEXTE spectral analysis 

Because of the non-imaging nature of the two main RXTE 
instruments the (total) source-flux estimation relies on accurate 
time-dependent instrumental and celestial background measure- 
ments. Although such models exist for the PCA, the complex- 
ity of the near environment of the AXPs in our sample makes 
it very difficult to derive reliable unbiased total flux estimates 
with the PCA in the 2-30 keV range. Specifically, all AXPs 
are located in a narrow strip along the Galactic plane where 
(large) gradients in the Galactic ridge emission exist (Valinia 
et al. 1998); IE 1841-045 and IE 2259 + 586 are located in 
supernova remnants; time- variable strong sources are present 
near to 4U 0142+61 (Be X-ray binary RX J0146.9+6121), IE 
1048. 1 - 5937 (the enigmatic r\ Carina) and 1RXS J 1708 - 4009 
(e.g. the strong and highly variable X-ray binaries OAO 1657- 
415 and 4U 1700-377). 

The rocking strategy applied during HEXTE operations in 
principle provides instantaneous background measurements, but 
also in this case the gradient in the Galactic ridge emission and 
the possible presence of other (strong) sources in both the on 
and off-source pointings (very serious for e.g. 1RXS J1708- 
4009 and 4U 0142 + 61) can result in unreliable (background 
subtracted) total-source-flux measurements. Therefore we aban- 
doned, in contrast to the work presented in Kuiper et al. (2004) 
for IE 1841-045 in Kes 73, the derivation of the total-source 
flux with the non-imaging PCA and HEXTE instruments. 

In this work we concentrate on the derivation of the time- 
averaged pulsed PCA/HEXTE spectra of the AXPs in our sam- 
ple. This can be done by determining the number of pulsed 
counts in differential PCA/HEXTE energy bands by fitting a 



truncated Fourier series 

N 

N(cf)) = oq + ^ at cos(2wk(j>) + bi sin(27r&<?!>) ( 1 ) 

k=l 

with (f> the pulse phase, to the measured pulse phase distribu- 
tions N((p). It turned out that 3 to 5 harmonics (N = 3/5) were 
sufficient to describe the measured distributions accurately for 
all energy intervals and AXPs in our sample. In the case of 
the PCA we derived for each PCU the energy response ma- 
trix (energy redistribution including the sensitive area) for the 
combination of observations listed in Table 1 and subsequently 
took the different PCU (screened) exposure times into account 
in the construction of the weighted PCU-combined energy re- 
sponse. The pulsed (excess) counts per energy band are fitted 
in a procedure assuming either an absorbed power-law (F 7 = 
K-e~ Nr, a -E~ T ), or an absorbed double power-law (F 7 = e~ Nli ' a ■ 
(K\ ■ E~ r ' +K2 ■ E~ Fl )) or an absorbed black body plus power- 
law (F 7 = e~ N "- a ■ (K ■ E~ r + Kbb ■ £ 7 /(exp(£ 7 /fer) - 1))) photon 
spectrum folded through the PCU-combined energy response. 

In the spectral "deconvolution" process of the HEXTE to- 
tal pulsed counts in almost all cases 4 the on-axis cluster A and 
B energy response matrices have been employed taking into ac- 
count the (slightly) different screened on-source exposure times 
for each cluster. The exposure times have been corrected for the 
considerable deadtime effects. 

3.3. INTEGRAL timing analysis 

The first step in an INTEGRAL timing analysis is to obtain 
a set of science windows for which the angular distance be- 
tween instrument pointing direction and target is within 14.°5 
to ensure that (a part of) the detector plane is illuminated by the 
target. The resulting list is further screened on erratic (ISGRI) 
count rate variations, indicative for particle effects due to Earth 
radiation belt passages or solar flare activities. These science 
windows are excluded for further analysis. Next, only events 
with rise times between 7 and 90 (see Lebrun et al. 2003, for 
definition), detected in non-noisy ISGRI detector pixels which 
have an illumination factor of more than 25% (i.e. at least 25% 
of a detector pixel must have been illuminated by the target) are 
passed for further analysis. The on-board event time stamps are 
corrected for known instrumental (fixed), ground station and 
general time delays in the on-board time vs. TT (Terrestrial 
Time) correlation (see e.g. Walter et al. 2003). The resulting 
event times in TT of the selected events are barycentered (using 
the JPL DE200 solar system ephemeris) adopting the Chandra 
X-ray positions of the AXPs and the instantaneous INTEGRAL 
orbit information. These barycentered events are finally folded 
using an appropriate timing model (v^v^v and the epoch) to 
yield pulse phase distributions for different energy bands be- 
tween 20 and 300 keV. The timing models (phase connecting 
ephemerides) are based on publicly available RXTE monitoring 
data of AXPs. Ephemerides have been generated for two AXPs 
in our sample, 1RXS J1708- 4009 and IE 1841-045, because 
at the INTEGRAL epoch (MJD > 52668 and MJD > 52698 
for 1RXS J1708-4009 and IE 1841-045, respectively) these 
were not available from existing literature (see Table 3; all with 
very small RMS values, required for extracting the weak pulsed 
signals). 

4 For HEXTE pointings with the target (AXP) far off-axis e.g. 4U 0142 + 61 
during an observation of HMXB RX J0146.9+6121, we took the reduction in 
the effective sensitive area due to the collimator response into account 
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Table 3 

Phase coherent ephemerides for 1RXS J1708 - 4009 and IE 1841-045, derived from RXTE PCA monitoring data and 

VALID FOR THE ANALYZED INTEGRAL OBSERVATIONS. 



AXP 


Start 


End 


Epoch 


V 


V 


V 


RMS 




[MJD] 


[MJD] 


[MJD,TDB] 


[Hz] 


X Kr 13 Hz/s 


X Kr 22 Hz/s 2 




1RXS J 1708 -4009 


52590 


52939 


52590.0 


0.09089812328(36) 


-1.59836(30) 


O.OO(Fixed) 


0.012 


IE 1841-045 


52726 


52982 


52726.0 


0.0848972590(50) 


-3.2217(85) 


4.22(83) 


0.015 



3.4. INTEGRAL spectral analysis 

The INTEGRAL IBIS ISGRI spectral analysis applied in our 
study is based on OS A 4.1 programs producing sky mosaics 
for the combination of many science windows in different en- 
ergy bands (Goldwurm et al. 2003). The resulting dead-time 
corrected ISGRI source-count rates per energy band are refer- 
enced to count rates measured from the Crab in similar energy 
bands (Crab calibration observations during INTEGRAL Revs. 
102/103 were used). 

We used for the total Crab photon emission spectrum the fol- 
lowing spectrum derived by Willingale et al. (2001) based on 
XMM-Newton observations of the Crab at energies between 0.3 
and 10 keV: photon index T = 2.108(6) with a normalization at 
1 keV of 9.59(5) ph/(cm 2 s keV). 

We verified the validity of the extrapolation of this spec- 
trum to energies between 15 and 250 keV using RXTE HEXTE 
Crab data. The HEXTE data utilized is from a long dedicated 
RXTE Crab observation (obs. id. 40805), performed between 
17-31 March 1999 and 18-19 Dec. 1999, yielding dead time 
corrected cluster and 1 exposures of 22.7 and 23.8 ks, re- 
spectively. Applying an overall (energy independent) normal- 
ization factor of 1.087 to be consistent with the BeppoSAX 
LECS/MECS/PDS spectrum (Kuiper et al. 2001) the derived 
15-250 keV HEXTE Crab total spectrum connects smoothly to 
the 0.3-10 keV Crab total spectrum derived by Willingale et al. 
(2001). The HEXTE Crab total spectrum can best be described 
by a power-law model with an energy dependent photon index 
as given below: 

F 7 = 1.5703(14) x (jB7 /o.06 3 55)- 2 097(2) - 0082(16) - ln(£ -/°- 06355) 

(2) 

In Eq.2, F 7 is expressed in ph/(cm 2 s MeV) and £ 7 in MeV. Up 
to ~ 180-200 keV the extrapolation of the spectrum given by 
Willingale et al. (2001) overlaps the HEXTE spectrum, and thus 
provides a proper representation for the ISGRI energy range we 
are effectively dealing with. 

Therefore, our method for deriving ISGRI spectra enables 
us to determine the 20-300 keV AXP spectra without detailed 
knowledge of the instrument energy response, which at the time 
of our analysis still contained significant uncertainties. Thus, 
we derive from these spatial analyses total source spectra (sum 
of pulsed and unpulsed components). 

For the construction of spectra for the pulsed components 
we extracted count rates from the phase distributions similarly 
as was done for the RXTE light curves. For the conversion 
to pulsed fluxes, also in this case, the Crab pulsed signal in 
the HEXTE 15 - 250 keV range (obs. id. 40805) was used as 
a reference source. This pulsed HEXTE spectrum is properly 
described by: 

F 7 = 0.4693(21) x (jB7 / .04844)- L955(7) - 0710(78) ln( ^/ 04844) 

(3) 



over the entire 15-250 keV energy range. This model has been 
verified on the pulsed ISGRI spectrum of PSR B1509-58 data, 
yielding consistent results as reported in Kuiper et al. (1999). 

3.5. COMPTEL spatial analysis 

During the long mission lifetime of CGRO (April, 1991 — 
May, 2000), most of the sky, particularly the Galactic Plane, 
has been viewed with long exposures. We used the exposure 
accumulated for each source over the total mission duration, 
amounting for 1RXS J1708-4009 5.3 Ms, 4U 0142 + 61 4.2 
Ms, IE 2259 + 586 4.9 Ms, IE 1048.1-5937 4.8 Ms, and for 
IE 1841-045 4.2 Ms. Skymaps and source parameters can be 
derived with the maximum likelihood method, which is imple- 
mented in the standard COMPTEL data analysis package. See 
for the implementation of the maximum likelihood method for 
data from a Compton telescope (de Boer et al. 1992), and for the 
specific treatment of the instrumental background structure in 
the COMPTEL data space (Bloemen et al. 1994). For the anal- 
ysis and data selections we followed the approach described by 
Zhang et al. (2004). Standard energy intervals were selected for 
the analysis: 0.75 - 3 MeV, 3-10 MeV and 10 - 30 MeV. For 
each of these energy intervals and for each AXP the source flux 
(or upper limit) was measured at the known source position in 
the maximum likelihood skymaps. 

4. 1RXS J1708-4009 

1 RXS J 1 708 - 4009 was discovered in 1 996 during the ROS AT 
(0.1-2.4 keV) all sky survey. X-ray pulsations at a 1 1-s period 
were subsequently detected with ASCA (Sugizaki et al. 1997). 
Its 0.8-10 keV spectrum turned out to be very soft and its gen- 
eral X-ray properties pointed to an AXP membership (see e.g. 
Israel et al. 1999a; Kaspi et al. 1999). From regular RXTE 
PCA monitoring observations performed between Jan. 1998 
and June 1999 a phase-coherent timing solution was obtained 
by Kaspi et al. (1999). This demonstrated a high level of rota- 
tion stability. Since these early monitoring observations, how- 
ever, the source experienced two glitches - one in Sept./Oct. 
1999 and a second in April 2001 - each with different recovery 
behaviour (Kaspi et al. 2000; Kaspi & Gavriil 2003; Dall'Osso 
et al. 2003). For the period between the two glitches Gavriil & 
Kaspi (2002) presented a phase-coherent timing solution with a 
positive i) indicative for a long-term glitch recovery. 

The morphology of the X-ray pulse profile of 1RXS J 1708 - 
4009 is changing as a function of energy (e.g. Sugizaki et al. 
1997; Israel et al. 2001; Gavriil & Kaspi 2002). Phase-resolved 
spectral analyses indeed showed significant spectral variations 
with pulse-phase, most pronounced in the photon power-law 
index (e.g. Israel et al. 2001; Rea et al. 2003, 2005). Further- 
more, the total phase-averaged unabsorbed 0.5-10 keV X-ray 
flux and photon spectral index appear to be time variable in a 
correlated way, with maximal fluxes and hardest spectra near 
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FIG. 1.— RXTE PCA/HEXTE pulse profiles of 1RXS J1708-4009 for 
energies in the range 1.8-222.9 keV combining data collected between 12 Jan 
1998 and 26 Oct 2003 (see Table 1). Two cycles are shown for clarity. The 
vertical dotted lines at phases 0.25 and 0.55 serve as a guide to the eye for 
alignment comparisons. Note the drastic morphology changes with energy. 



the two glitch epochs (Rea et al. 2005). 

At optical/IR wavelengths 2 potential counterparts were iden- 
tified within the Chandra O'.'l HRC-I error circle (see e.g. Israel 
et al. 2003; Safi-Harb & West 2005, for more details). A search 
for radio emission at 1.4 GHz from 1RXS J 1708 -4009 only 
yielded a 5a upper limit of 3 mJy at the position of the AXP 
(Gaensler et al. 2001). 

Given the softness of the 0.5-10 keV X-ray spectra, the IN- 
TEGRAL detection reported by Revnivtsev et al. (2004) of a 
point source at the position of 1RXS J 1708 -4009 between 18 
and 60 keV was a big surprise. Below we will present in detail 
the new high-energy characteristics of this AXP derived in this 
work: a) the discovery of the pulsed emission above ~ 10 keV 
(profiles, spectra) using RXTE PCA/HEXTE and IBIS ISGRI 
data; b) ISGRI and COMPTEL results on the total emission. 

4.1. 1RXS Jl 708 - 4009 timing characteristics 

4.1.1. RXTE PCA/HEXTE pulse profiles 

Applying the timing analysis procedures outlined in Sect. 3.1 
to the full set of RXTE observations of 1RXS J 1708 -4009 
listed in Table 1, resulted in a compilation of high-statistics 
time-averaged PCA/HEXTE pulse profiles for energies between 
~ 2-220 keV (see Fig. 1). The ephemerides used in the fold- 
ing/correlation process (see Sect. 3.1) are given in Kaspi et al. 
(2000); Gavriil & Kaspi (2002); Kaspi & Gavriil (2003). For 
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FIG. 2. — Phase angles as a function of energy for the first 3 harmonics used 
in the truncated Fourier series fit (see Eq. 1) of the RXTE PCA pulse profiles 
of 1RXS J1708-4009. The harmonics are labeled with their corresponding 
number. 



the first time pulsed emission is detected above ^10 keV: the 
non-uniformity significance of the 16.1-32.0 keV PCA pulse 
phase distribution (see Fig. Id) is 14.2cr applying a Zj-test and 
the HEXTE 35.2-222.9 keV profile (see Fig. If) deviates from 
uniformity at a 5.2tr level. Above 35.2 keV the significances 
in the HEXTE 35.2-64.1 and 74.3- 222.9 keV bands (the inter- 
mediate energy window with a large instrumental background 
feature has been omitted) are both 3.75c. Drastic morphol- 
ogy changes with energy are visible. The decomposition of the 
pulse profiles in terms of a finite number of harmonics (see Eq. 
1) provides a means to visualize a change in morphology with 
energy. The power {a 2 k + b 2 k , see Eq.l), derived from the time- 
averaged PCA pulse profiles of 1RXS J 1708 -4009, in the first 
harmonic is dominant over the power in the second and third 
harmonics, and the power in harmonics with k > 4 can be ne- 
glected (see also Gavriil & Kaspi 2002). From Eq. 1 one can 
define a phase angle = arctan^/^) for each harmonic k. 
The energy dependence of for the first three harmonics is 
shown in Fig. 2. For each harmonic it reveals a very smooth 
variation of the phase angle with energy. The shape of the pro- 
file is changing drastically between 2 and 10 keV. For energies 
above ^15 keV the phase angles for the 3 considered harmon- 
ics seem to converge to constant values, a necessary condition 
for stable pulse shapes. 

4. 1 .2. INTEGRAL IBIS ISGRI pulse profiles 

We also performed a timing analysis for 1RXS J 1708 -4009 
using IBIS ISGRI data. Data from science windows taken dur- 
ing INTEGRAL Revs. 36-120 satisfying our 14.°5 off-axis con- 
straint were included (effective on-axis exposure after screen- 
ing <~ 1,360 ks). The processing followed the guidelines pre- 
sented in Sect. 3.3 using the 1RXS J1708-4009 ephemeris 
generated from RXTE monitoring observations, given in Ta- 
ble 3. In the integral 20-300 keV ISGRI band we obtained a 
non-uniformity significance of 5.9cr applying a Z| test, which 
is comparable to the HEXTE result. In differential energy bands 
we found: 20-75 keV, 4.3cr, and 75-300 keV 3.6a (see Fig. 3 
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IBIS ISGRI: 20-75 kcV a 




IBIS ISGRI: 75-300 kcV b 




0.0 0.5 1.0 1.5 2.0 

Phase 

FIG. 3.— INTEGRAL IBIS ISGRI pulse profiles of 1RXS J1708 -4009 for 
two energy ranges. The non-uniformity significances are 4.3<r and 3.6a for 
20-75 keV and 75-300 keV, respectively. Pulse maxima are found near phase 
~ 0.2, corresponding to phase 0.55 in Fig. 1. 

for the corresponding pulse profiles). The HEXTE and the IS- 
GRI profiles above 75 keV are very similar, and suggest that 
the hard-X-ray 1RXS J 1708 -4009 profile exhibits less struc- 
ture than found below 10 keV. From these initial ISGRI timing 
results it is clear that highly significant profiles can be expected 
in the near future when significantly more IBIS ISGRI data on 
this source become available. 

4.2. 1RXS Jl 708 - 4009 spectral characteristics 

In this section we present new high-energy spectral informa- 
tion above 2.5 keV up to 30 MeV for 1RXS J1708 - 4009: (a) 
(time-averaged) pulsed emission from RXTE PCA and HEXTE; 
(b) pulsed emission from INTEGRAL IBIS ISGRI; (c) total 
(pulsed and unpulsed) emission from ISGRI and upper lim- 
its to the total emission from CGRO COMPTEL. Finally, the 
new spectra are compared with spectra reported earlier for en- 
ergies below 10 keV (Rea et al. 2003, 2005, for BeppoSAX 
LECS/MECS 0.4-10.8 keV and XMM Newton MOS/PN 0.5- 
10 keV, respectively). 

4.2.1. RXTE PCA/HEXTE pulsed spectrum 

The spectral procedures employed for the RXTE PCA and 
HEXTE data (see Sect. 3.2) resulted in a high-statistics deter- 
mination of the spectrum of the time-averaged pulsed emission 
of 1RXS J1708-4009 in the - 2.5-220 keV energy range. 
The PCA (aqua) flux values are derived assuming an absorbed 
double power-law spectral model, and are shown in a vF v repre- 
sentation in Fig. 4. Also drawn is the best fitting spectral model 
to the PCA data points (2.5-36.9 keV; ^=1.11 for 12 de- 
grees of freedom; dashed line). The assumed absorbing Hydro- 
gen column density Nh in the spectral fit was 1.36 x 10 22 cm~ 2 
(Rea et al. 2003). The two power-law components become 
equally strong at E cross = 21.7 ±2.4 keV: below this energy 
the power-law component with index Y\ = 2.60 ±0.01 domi- 
nates and above a component with very hard spectrum, index 
T2 = -0.12 ±0.07. It is clear that the pulsed spectrum hardens 
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FIG. 4. — A vF v spectral representation of the total and pulsed high-energy 
emission from 1RXS J1708-4009. The aqua (PCA), blue (HEXTE) and ma- 
genta (IBIS ISGRI) data points show the time-averaged 2.5-300 keV pulsed 
spectrum. The black dashed line shows the best power-law model fit to the 
combination of PCA, HEXTE and ISGRI pulsed flux values for energies above 
~ 15 keV. The other measurements refer to the total emission spectrum: 0.5-10 
keV, BeppoSAX LECS/MECS and XMM Newton spectral models at different 
epochs (Rea et al. 2003,2005); 20-300 keV time-averaged (Revs. 36-106) 
IBIS ISGRI spectrum; and 0.75-30 MeV, time-averaged CGRO COMPTEL 
2cr upper-limits. Note the drastic hardening of the pulsed spectrum near 20 
keV. The COMPTEL upper limits require another spectral break somewhere 
between 300 and 750 keV. 

dramatically above 20 keV, however, the spectrum has to soften 
considerably to be consistent with the index of 1 .01(12) derived 
from a combination of PCA, HEXTE and IBIS ISGRI pulsed 
flux measurements for energies > 15 keV (see Sect.4.2.2). The 
HEXTE flux values (simultaneously derived) are fully consis- 
tent with the upturn found in the PCA spectrum (see Fig.4, blue 
datapoints). 

4.2.2. INTEGRAL IBIS ISGRI pulsed spectrum 

Pulse profiles for 1RXS J 1708 - 4009 have been generated in 
three energy intervals from 20 to 300 keV using IBIS ISGRI 
data from Revs 36-120. The pulsed (excess) count rates were 
determined and transformed to pulsed flux values, calibrated 
on the pulsed Crab spectrum. The three flux values are also 
included in Fig. 4. They are consistent with the data points 
derived from RXTE HEXTE for the same energy window, but 
measured at different epochs. Within the present statistical ac- 
curacies there is no indication for long-term flux variability 
above 20 keV, contrary to the reported strong flux variability 
for energies below 10 keV, clearly visible in Fig. 4. A power- 
law fit through the PCA, HEXTE and IBIS ISGRI pulsed flux 
measurements for energies above 15 keV yielded a photon in- 
dex of 1.01 ±0.12. This model is shown in Fig. 4 as a dashed 
black line for the 15-300 keV energy range. There is no indica- 
tion yet for a second spectral break up to 300 keV. 

4.2.3. INTEGRAL IBIS ISGRI and COMPTEL total spectrum 

1RXS J1708-4009 was detected at a 6.5cr level in the 18-60 
keV energy band by Revnivtsev et al. (2004) using IBIS ISGRI 
data from a 2 Ms ultra deep INTEGRAL survey of the Galac- 
tic center region performed in Aug. - Sept. 2003 (INTEGRAL 
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Revs. 104-107 & 111-113). Its 18-60 keV flux was 2.2 ±0.3 
mCrab. In this work we analyzed (almost independent) IBIS IS- 
GRI data from all publicly available INTEGRAL observations 
performed between revolutions 36 and 106 in which the angu- 
lar distance between 1RXS J 1708 -4009 and the science win- 
dow pointing was < 14.°5. The effective on-axis exposure after 
screening is about 974 ks. Mosaic images have been made us- 
ing OS A 4. 1 procedures for 5 broad differential energy bands, 
20-35, 35-60, 60-100, 100-175 and 175-300 keV, in order to 
determine the total hard X-ray/soft 7-ray spectrum of this AXR 
1RXS J1708-4009 was detected significantly only in the 20- 
35 and 35-60 keV energy bands with flux levels consistent with 
the 18-60 keV flux measurement by Revnivtsev et al. (2004). 
The resulting IBIS ISGRI spectral measurements are shown as 
purple data points in Fig. 4. Unfortunately, the statistics in 
these IBIS ISGRI data points are still too poor to constrain the 
underlying emission model in the 20-300 keV range: the fit 
quality specified by a reduced of 1.77 for 3 d.o.f. assum- 
ing a simple power-law model is still acceptable; the resulting 
photon index is 1.44 ±0.45, slightly softer than, but consistent 
with, the pulsed flux spectrum derived from the combined fit to 
the PCA, HEXTE and IBIS ISGRI pulsed flux measurements 
at > 15 keV. 

Finally, at energies above 750 keV till 30 MeV we generated 
skymaps using all CGRO COMPTEL observations with 1RXS 
J 1708 -4009 in its field of view, spread over its full 9 year mis- 
sion lifetime. The source was not seen for any of the standard 
energy intervals, and we derived (2a) flux upper limits, also 
shown in Fig. 4. 

4.2.4. Discussion of spectra and pulsed fraction oflRXS 
J 1708 -4009 

Comparing the total IBIS ISGRI spectrum derived in the spa- 
tial analysis with the pulsed time-averaged spectra from RXTE 
PCA/HEXTE and ISGRI, we see that the 1RXS J 1708 -4009 
emission is consistent with about 100% pulsation for energies 
in excess of 20 keV. For energies below 20 keV, the situation is 
very different. To highlight this, we included in Fig.4 the emis- 
sion models (best spectral fits) for the total 0.5-10 keV X-ray 
spectra measured at two different epochs (Rea et al. 2003, 2005, 
for BeppoSAX LECS/MECS 0.4-10.8 keV and XMM Newton 
MOS/PN 0.5-10 keV, respectively). Given the apparent drastic 
spectral change in both shape and normalization for energies 
below 10 keV we cannot uniquely quantify the pulsed fraction 
as a function of energy comparing with the time-averaged high- 
statistics PCA pulsed flux measurements. Obviously, the Bep- 
poSAX measurement during August 2001 does not connect to 
the RXTE/INTEGRAL spectra. 1RXS J 1708 -4009 was ap- 
parently in a (very) different high state during the BeppoSAX 
observation, which took place in the recovering phase from the 
secondly reported glitch (see Rea et al. 2005). The XMM New- 
ton spectrum connects very smoothly to the total spectrum mea- 
sured by INTEGRAL, and suggests a variation in pulsed frac- 
tion from ~ 25% at 2 keV to 100% at 20 keV. Our time av- 
eraged results for the total and pulsed emission above 20 keV 
suggest a very stable behaviour of this AXP at these higher en- 
ergies. It seems that a stable and very hard pulsed component 
is dominating the emission at these energies. Interestingly, the 
COMPTEL MeV upper limits then require a spectral break in 
this hard spectrum of 1RXS J 1708 -4009 somewhere between 
300 and 750 keV. 



5. 4U 0142 + 61 

4U 0142 + 61 was discovered in the early seventies by the 
scanning Uhuru X-ray observatory (see e.g. Forman et al. 1978). 
The X-ray sky survey performed by the SSI (1.5-20 keV) on 
Ariel V (Warwick et al. 198 1) and X-ray observations by SAS-3 
and HEAO-1 confirmed the source and refined its position con- 
siderably (see e.g. Reid et al. 1980, for an X-ray position with 
23" accuracy). Data from the non-imaging ME collimator (45' 
FWHM) aboard EXOS AT revealed a coherent 25 minute signal, 
which was later absent in two shorter EXOSAT observations 
(White et al. 1987). These authors also presented more accu- 
rate (3"2) positional information for 4U 0142+61 using both 
EXOSAT LEIT and Einstein (HEAO-2) HRI data and demon- 
strated the lack of a bright optical counterpart, ruling out the 
presence of a massive companion. 

The discovery in ROSAT all-sky survey data of a Be/X-ray 
binary, RX J0146.9+6121, by Motch et al. (1991) at only 24' 
from 4U 0142 + 61 showed that there had been source confu- 
sion. Revisiting the EXOSAT 1985 data, Israel et al. (1994) 
discovered a 8.7 s periodicity, only detectable in the 1-3 keV en- 
ergy band. Finally, using ROSAT PSPC data, Hellier (1994) un- 
ambiguously demonstrated that the 8.7 s periodicity originated 
from 4U 0142 + 61 and the 25 min signal from the Be/X-ray bi- 
nary RX J0146. 9+6121, which moreover showed transient ac- 
tivity. 

A good quality phase-averaged 0.5-10 keV X-ray spectrum 
of 4U 0142 + 61 was derived by White et al. (1996) analyz- 
ing ASCA SIS/GIS data. The physically most pausible model 
consists of an absorbed black-body plus power-law component 
with best-fit model parameters: black-body temperature 0.386(5) 
keV; photon-index 3.67(9) and N H 9.5(4) x 10 21 cm" 2 . These 
authors, using data covering a much wider time baseline, also 
confirmed the spin-down timescale of <~ 1.2 x 10 5 year, re- 
ported earlier by Hellier (1994). The spin-down energy release 
is much lower than the observed X-ray luminosity, thereby ex- 
cluding 4U 0142 + 61 to be a spin-down powered pulsar. 

RXTE monitoring observations since November 1996 (Gavriil 
& Kaspi 2002) proved that 4U 0142 + 61 is a very stable ro- 
tator, like spin-down powered pulsars, and ASCA data taken 
in a 2 year gap of RXTE observations seemed to indicate a 
pulsar-like glitch in its rotation behaviour (Morii et al. 2005). 
Time-averaged RXTE PCA pulse profiles in two different en- 
ergy bands showed significant morphology changes with en- 
ergy, consistent with earlier results of much lower significance 
(Israel et al. 1994; White et al. 1996; Israel et al. 1999b). 

Recent X-ray observations of 4U 0142 + 61 by Chandra (Juett 
et al. 2002; Patel et al. 2003) and XMM (Gohler et al. 2005) im- 
proved not only the positional accuracy down to 0."5, but also 
provided high-quality (phase-resolved) spectra which are mu- 
tually consistent and appeared featureless. The phase-resolved 
spectra showed significant changes in photon-index and black 
body temperature as a function of pulse phase. 

At lower energies, Hulleman et al. (2000) discovered a faint 
optical counterpart within the 4U 0142 + 61 Einstein X-ray er- 
ror circle, followed by the detection of optical pulsations with 
a high pulsed fraction from this counterpart (Kern & Martin 
2002). No radio counterpart has been found yet (Gaensler et al. 
2001). 

At higher energies the discovery was reported by den Har- 
tog et al. (2004) of hard X-rays/soft 7-rays (20 - 150 keV) in 
INTEGRAL IBIS ISGRI data of a deep 1.6 Ms observation of 
the Cassiopeia region. A detailed analysis of this long obser- 
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FIG. 5.— ASCA GIS, RXTE PCA and HEXTE pulse profiles of 4U 
0142 + 61 for energies in the range 0.5-156.8 keV. Two cycles are shown for 
clarity. The RXTE profiles shown in panels d-i are time-averaged and based on 
observations performed between March 1996 and September 2003. The GIS 
profiles (panels a-c) are from ASCA observations performed in July/August 
1999 totaling ~ 120 ks of exposure. For the first time significant deviations 
from uniformity are visible at energies > 8 keV: Z| is 16.6<r and 5.9<r for 
the PCA 8.0-16.3 and 16.3-31.8 keV bands, respectively; the significances for 
the HEXTE profiles are 3.4<r and 2.0a for the 21.2-50.3 and 50.3-156.6 keV 
bands. Note the drastic morphology changes with energy. 



vation (den Hartog et al. 2006), proved that the spectrum in 
this energy range is very hard with power-law photon index 
T = 0.73 ± 0.17. The X-ray luminosity between 20 and 100 
keV was found to be 5.9 x 10 34 erg s 1 , a factor 440 higher than 
the rotational energy loss. They also reported flux upper limits 
from COMPTEL between 0.75 and 30 MeV, which indicated 
that the hard spectrum has to break/bend between ^100 keV 
and 0.75 MeV. 

In this work we will present for 4U 0142 + 61: a) for the first 
time the PCA/ HEXTE timing and spectral results for energies 
above ~ 8 keV using all publicly available RXTE data, b) tim- 
ing and spectral results from analysis of ASCA GIS 0.5 - 10 
keV data, and c) a spectrum with higher statistical accuracy of 
the total emission as seen by IBIS ISGRI over the 20-300 keV 
energy range using all available public, Core Program (Galac- 
tic plane scans; 4U 01 15+63 Target of Oppotunity Observation) 
and Open Time (Cassiopeia region; IGR J0029 1+5934 TOO) 
INTEGRAL data. 

5.1. 4U 0142 + 61 timing characteristics 

Time-averaged RXTE pulse profiles of 4U 0142 + 61 are shown 
in Fig. 5 for energies between 2.2 and 156.6 keV (PCA/HEXTE). 
We used in the folding/correlation process (see Sect. 3.1) the 
ephemerides given by Gavriil & Kaspi (2002). The 2.2-4.0 keV 
PCA profile (panel d) shows two distinct pulses near phases 0.2 
and 0.55. Moving up in energy the pulse near 0.2 loses sig- 
nificance and is gone for energies above ~ 8 keV (panels e-f). 




FIG . 6 . — Phase angles as a function of energy for the first 3 harmonics of the 
truncated Fourier series fit of the RXTE PCA pulse profiles of 4U 0142 + 61. 
Note the sign reversal of the energy derivative of & a for the second harmonic 
near 5 keV. 



Instead above ~ 8 keV a feature near phase 0.9 pops up, which 
is visible up to <~ 50 keV (panels f-h). Above 50 keV (panel i) 
the phase region between the two pulses at 0.55 and 0.9 seems 
to be filled in by a new component, but the statistics is too poor 
to make stringent conclusions. 

Most importantly, we detect for the first time significant pulsed 
emission from4U 0142 + 61 at energies above ~ 8 keV: 16.6<r 
and 5.9a for the PCA 8.0-16.3 and 16.3-31.8 keV bands, re- 
spectively. To investigate further the increase in strength of 
the pulse near phase 0.2 towards lower energies, we extended 
the energy baseline by including also our results from a timing 
analysis of ASCA GIS (0.5- 10 keV) data from observations per- 
formed in July/August 1999 (~ 120 ks exposure). Three pro- 
files (energy bands 0.5-1.7, 1.7-3.0 and 3.0-10 keV) are shown 
in panels a-c of Fig. 5. Indeed, for the pulse near phase 0.2 the 
trend of increasing strength towards lower energies is continued 
in the ASCA GIS 0.5-1.7 and 1.7-3.0 keV energy ranges. 

A decomposition of the PCA pulse profiles in terms of 5 
Fourier harmonics shows that above 2 keV most of the signal 
power is embedded in the first harmonic which becomes more 
and more dominant over the power in the second harmonic, the 
next powerful harmonic, with increasing energy. Below ~ 2 
keV the power in the second harmonic is dominant over that in 
the first. The behaviour of the phase angles as a function of 
energy for the first three harmonics is shown in Fig. 6. Simi- 
lar to 1RXS J 1708 -4009 but with different trends, also now a 
smooth energy dependency of is shown for each harmonic 
with significant power. The exception is the derivative of <fr 2 a to 
energy which flips sign near 5 keV. For 4U 0142 + 61 it is less 
obvious that the pulse profile remains stable above 20 keV, than 
seems to be the case for 1RXS J1708 - 4009. 

5.2. 4U 0142 + 61 spectral characteristics 

In this section we present: a) the time-averaged pulsed spec- 
trum (2.2-102 keV) of4U 0142 + 61 based on RXTE PCA and 
HEXTE measurements. The pulsed spectrum is extended down 
to 0.8 keV by including our spectral results from the ASCA 
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FIG. 7. — A vF v spectral representation of the total and pulsed high-energy 
emission from 4U 0142 + 61. The aqua, blue and dark cyan data points/curves 
represent the pulsed emission component of the spectrum (0.8-102 keV) based 
on RXTE PCA/HEXTE observations (time-averaged) and ASCA GIS. Note 
the dramatic hardening of the pulsed spectrum near 10 keV. The other mea- 
surements refer to the total emission spectrum: 0.8-10 keV, Chandra ACIS 
(Patel et al. 2003), 20-300 keV, INTEGRAL IBIS ISGRI (this work) and 0.75- 
30 MeV, CGRO COMPTEL (den Hartog et al. 2006). The COMPTEL 2a 
upper limits require a spectral break somewhere between 140 and 750 keV. 



GIS observations (~120 ks); b) the time-averaged total (pulsed 
and unpulsed) spectrum from our analysis of INTEGRAL IBIS 
ISGRI skymaps (20-300 keV), in comparison with that from 
Chandra ACIS-S CC-mode data (0.5-7 keV; Patel et al. 2003) 
and the CGRO COMPTEL (0.75-30 MeV) upper limits (den 
Hartog et al. 2006). 

5.2.1. RXTE PCA/HEXTE and ASCA GIS pulsed spectrum 

Assuming an underlying absorbed {Nh = 9.3 x 10 21 cm" 2 ; Pa- 
tel et al. 2003) double power-law photon model we fitted the 
pulsed PCA excess counts (2.2-31.5 keV) in a forward folding 
procedure (see Sect. 3.2), resulting in a high-statistics time- 
averaged pulsed spectrum. The best model fit (x 2 . = 0.7; d.o.f. 
15 - 4) yielded for the soft and hard photon indices Ti = 4.09 ± 
0.04 and T2 = -0.8^ 07, respectively. This is shown as a dashed 
aqua curve in Fig. 7, together with the "deconvolved" unab- 
sorbed flux values (aqua colored open circles). The power-law 
model components become equally strong at 1 1 .46 ± 0.64 keV. 

The striking spectral hardening of the pulsed emission of 4U 
0142 + 61 around 10 keV is dramatic, much more pronounced 
than the hardening observed for 1RXS J1708 -4009. 

The HEXTE pulsed flux values (Fig. 7; dark-blue open squares) 
are in line with the extrapolation of the PCA hard-power-law 
model component. 

The ASCA GIS 0.8-10 keV pulsed flux values which we de- 
rived from the July/ August 1999 observations are also shown in 
Fig. 7 (filled squares) and are within the systematic/statistical 
uncertainties consistent with the time-averaged PCA fluxes in 
the overlapping energy range. This suggests a rather stable 
pulsed spectrum both in shape and normalization. 

5.2.2. INTEGRAL IBIS ISGRI total spectrum 



den Hartog et al. (2006) report the detection and total spec- 
trum of 4U 0142 +-61 in the 20-100 keV energy range analyz- 
ing IBIS ISGRI data from INTEGRAL observations of the Cas- 
siopeia region performed in Dec. 2003 (Revs. 142 - 148). In 
this work all available INTEGRAL data (open time/core pro- 
gram) have been used from observations made between 3 March 
2003 (Rev. 47) and 28 December 2004 (Rev. 269) in which 
the source was within 14.°5 from the pointing axis. The total 
screened exposure time for the set of accumulated science win- 
dows is about 2 Ms, which corresponds to an effective on-axis 
exposure of about 858 ks. Image mosaics have been gener- 
ated using OS A 4. 1 software in 7 differential energy bands, 20- 
30, 30-45, 45-65, 65-95, 95-140, 140-205 and 205-300 keV. In 
these (time-averaged) maps 5 4U 0142 + 61 is clearly detected 
up to and including the 95-140 keV energy range (7.8 ±1.0 
mCrab). The 20-300 keV flux measurements are shown as pur- 
ple data points in Fig. 7. A simple power-law fit to these flux 
points yielded a hard photon index of 1.05 ±0.11 (x 2 = 0.86 
for d.o.f. 6-2; the purple dashed line in Fig. 7). The differ- 
ence with the value (0.73 ±0.17) reported for the index by den 
Hartog et al. (2006) is 1.6<r. 

5.2.3. Discussion of spectra and pulsed fractions of 4U 
0142 + 61 

The INTEGRAL spectrum of the total emission above 20 
keV can be compared with the same Chandra spectrum for en- 
ergies below 10 keV (see Fig. 7 the dark orange solid line; 
Patel et al. 2003). As is the case for the pulsed spectra above 
and below ^10 keV, these two spectra are drastically differ- 
ent, and together reveal a sharp minimum in luminosity of 4U 
0142 +-61 around 10 keV. 

Including now in this high-energy spectral picture the time- 
averaged CGRO COMPTEL 2a flux upper limits (red trian- 
gles in Fig. 7) from den Hartog et al. (2006) confirms that the 
power-law model statisfactorily describing the 20-300 keV to- 
tal spectrum will not extend into the MeV range, but must break 
somewhere between 140 and 750 keV. To investigate this fur- 
ther, we also fitted a spectral model with an energy-dependent 



photon index, F 1 = K E 



-(r+Q ln(£ T )) 



to the 20-300 keV IBIS IS- 



GRI measurements (see Fig. 7 dashed dotted purple line). The 
model extrapolation towards the MeV energy range is consis- 
tent with the COMPTEL upper limits, but the improvement of 
the fit is insufficient to claim a change in spectral shape in the 
20-300 keV window. 

Comparing now the total and pulsed high-energy spectra of 
4U 0142 +-61, the Chandra total flux values below 10 keV are 
about 10 times higher than the time-averaged RXTE PCA and 
ASCA GIS pulsed flux measurements. This is consistent with 
pulsed fraction estimates of ~ 10% reported by Patel et al. 
(2003) and Gohler et al. (2005) using solely Chandra and XMM 
Newton data, respectively. The total spectrum measured by IN- 
TEGRAL above 20 keV with slope 1.05 ±0.11 is significantly 
softer than the pulsed spectrum measured by RXTE above ~ 
10 keV with slope -0.8±[}q7. As a result, the pulsed fraction 



appears to vary with energy from > 
between 80 and 100 keV. 



10% at 20keVto- 100% 



6. IE 2259 + 586 

Near the center of SNR CTB 1 09 (G 1 09 . 1 - 1 .0) a strong com- 
pact X-ray source was found by Gregory & Fahlman (1980), 

5 Be X-ray binary RX J0146.9+6121 at only 24' from 4U 0142 + 61 was de- 
tected, fully resolved, only in the first three energy bands. 
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later dubbed IE 2259 + 586, which appeared to be an X-ray 
pulsar with a period of 3.489(2) s (Fahlman & Gregory 1981). 
Analyzing more IPC data, this period turned out to coincide 
with the first harmonic of a fundamental 6.978 s period (Fahlman 
& Gregory 1983). 

No radio enhancement at the position of the X-ray pulsar 
could be identified down to a level of 0.5 mJy (WSRT, 21 cm; 
Hughes et al. 1984), 0.2 mJy (VLA, 20 cm; Gregory et al. 1983) 
and even 50^Jy (VLA, 20cm; Coe et al. 1994; Kaspi et al. 
2003). On the other hand, the lack of a bright visual counterpart 
(V > 21) ruled out a supergiant or massive main-sequence star 
association (Fahlman & Gregory 1981). 

Further X-ray observations of IE 2259 + 586 in the eighties 
with Tenma (Astro-B), EXOSAT and Ginga (Astro-C) (see e.g. 
Koyama et al. 1987; Morini et al. 1988; Hanson et al. 1988; 
Koyama et al. 1989; Iwasawa et al. 1992) revealed a steady 
spin-down, too slow to power the observed X-ray luminosity, 
and a very soft spectrum. X-ray flux measurements from two 
different Ginga observations spread ~ 8 months apart indicated 
considerable flux variations, accompanied with clear pulse mor- 
phology changes and probably a decrease in spin-down rate 
(Iwasawa et al. 1992). The X-ray picture was refined consider- 
ably in the nineties using data from BBXRT, ASCA (Astro-D), 
ROSAT, BeppoSAX and RXTE. Corbet et al. (1995) showed 
for the first time that a two-component spectral model - black- 
body plus power-law - could fit the ASCA and BBXRT spec- 
tral data satisfactorily without invoking spectral (line) features. 
This finding was confirmed by Rho & Petre (1997) analyzing 
ROSAT PSPC, ASCA and BBXRT data simultaneously, and by 
Parmar et al. (1998) analyzing BeppoSAX LECS and MECS 
data. RXTE timing measurements (Mereghetti et al. 1998) fur- 
ther tied down earlier limits on the projected semi-axis to 0.03 
Its, leaving room only for a white dwarf companion, helium- 
burning star with mass smaller than 0.8 M or a main-sequence 
star viewed under very small inclination angles. 

Kaspi et al. (1999) obtained for the first time a phase co- 
herent timing solution for IE 2259 + 586, indicating great sta- 
bility (rms 0.01 cycles) over the 2.6 year timespan (29-Sept- 
1996 - 12-May-1999) of their RXTE monitoring observations. 
This stability makes binary accretion scenarios very unlikely 
and favours a magnetar interpretation. Additional RXTE moni- 
toring data showed that this rotational stability was maintained 
throughout an extended 4.5 yr period, although the inclusion of 
v was required in the timing model (Gavriil & Kaspi 2002). 
These authors also found that the pulse morphology did not 
change significantly with time (cf. Iwasawa et al. 1992, who 
claimed changes in morphology with time). Moreover, there 
was no evidence for large variability in the pulsed flux in line 
with earlier work by Baykal et al. (2000) using RXTE data, but 
in contrast with the Ginga findings (Iwasawa et al. 1992). Fur- 
thermore, Gavriil & Kaspi (2002) showed, convincingly for the 
first time, that the pulse profile morphology of IE 2259 + 586 
changes with energy (cf. Hanson et al. 1988, for earlier indica- 
tions). 

Observations with the Chandra X-ray observatory in January 
2000 provided an X-ray position with subarcsecond accuracy 
(Hulleman et al. 2001; Patel et al. 2001). In the 0"6 error circle 
(99%) a very faint (K s = 21 .7 ± 0.2 mag) near-infrared counter- 
part was identified (Hulleman et al. 2001), excluding models in 
which the source is powered by disk accretion. The total X-ray 
(0.5-7.0 keV) spectrum of IE 2259 + 586 appeared featureless 
and could best be described by a combination of a black body 



kT = 0.412(6) keV plus power-law T = 3.6(1) absorbed by a Hy- 
drogen column Nh = (9.3 ±0.3) x 10 21 cm" 2 (Patel et al. 2001). 

In June 2002 RXTE observed an outburst of IE 2259 + 586 
during which both the pulsed and persistent X-ray emission in- 
creased by more than an order of magnitude relative to their 
quiescent levels (Kaspi et al. 2003). During the course of the 
observation a spectral softening occured, and a significant pulse 
profile change was observed. In the meantime the pulsar under- 
went a sudden spin-up (glitch) followed by a large increase in 
spin-down rate lasting for more than 18 days. Furthermore, 80 
X-ray bursts were detected during the 14.4 ks RXTE observa- 
tion with durations ranging from 2 ms to 3 s. The outburst prop- 
erties of IE 2259 + 586 share strong similarities with SGR out- 
burst characteristics unifying thereby conclusively AXPs and 
SGRs and supporting strongly a magnetar interpretation. More 
detailed information on the burst characteristics can be found 
in Gavriil et al. (2004). Woods et al. (2004) presented a com- 
parison of the X-ray emission characteristics of IE 2259 + 586 
before, during and after the June 2002 outburst using data from 
XMM-Newton and RXTE. They quantified the changes of the 
temporal and spectral properties and derived recovery timescales 
The X-ray flux increase and subsequent decay can be described 
by two distinct components: one component linked to the burst 
activity with a timescale of ~ 2 days and a second component 
which decays over the course of the year according to a power- 
law in time (F oc t a ) with index a = -0.22 ±0.01. The lat- 
ter component behaves similarly in time as the observed near- 
infrared flux decay (Tam et al. 2004) and thus these seem to 
be linked to a common physical mechanism likely acting in the 
pulsar's magnetosphere. 

In this work we derived for the first time the time-averaged 
timing/spectral properties of IE 2259 + 586 for energies above 
~ 8 keV using archival RXTE PCA/HEXTE data and we com- 
pared the pulsed spectrum with the upper limits from IBIS IS- 
GRI and COMPTEL (den Hartog et al. 2006). 

6.1. IE 2259 + 586 timing characteristics 

The time-averaged RXTE PCA/HEXTE pulse profiles com- 
bining all observations of IE 2259 + 586 listed in Table 1, ex- 
cept the 14.4 ks period of outburst in 2002, are shown in Fig. 8 
for energies between 2.2 and 27 keV. In the folding/correlation 
process (see Sect. 3.1) we used the ephemeris given by Gavriil 
& Kaspi (2002). The screened PCU-2 exposure amounts ~ 
747 ks, while the corresponding dead-time corrected on-source 
HEXTE exposures are 256.4 and 267.7 ks for cluster-0 and 1, 
respectively. For the first time significant deviations from uni- 
formity are visible at energies > 8 keV: Z| is 10.6cr, 5.2c and 
3.1cr for the PCA 8.3-11.9, 11.9-16.3 and 16.3-24 keV bands, 
respectively (Fig. 8 panels c-e). The HEXTE sensitivity is still 
too low to detect the underlying pulsed emission above ^15 
keV (Fig. 8 panel f). 

The morphology of the double-peaked pulse profile changes 
gradually moving up in energy from 2.2 to 24 keV: the domi- 
nant pulse near phase 0.7 in Fig. 8 at energies below ~ 4 keV 
loses significance with respect to the second pulse near phase 
0.25. The latter dominates at energies above ~ 8 keV. Also note 
the existence of a narrow pulse-like feature in the deepest min- 
imum near phase 0.45, most striking in panel b (4-8.3 keV) of 
Fig. 8. The latter feature is also clearly visible in Fig. 7c of 
Gavriil & Kaspi (2002). 

More quantitative information on the energy dependency of 
the pulse morphology can be obtained through a decomposi- 
tion of the pulse profiles in Fourier components. At least 5 har- 
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Arbitrary Phase 

FIG. 8.— RXTE PCA/HEXTE pulse profile collage of IE 2259 + 586 for 
energies in the range 2.2-27.0 keV combining data collected between 29 Sept 
1996 and 28 Oct 2003 (see Table 1). Two cycles are shown for clarity. The 
vertical dotted lines at phases 0.25 and 0.70 serve as a guide to the eye for 
alignment comparisons. Clear pulse profile morphology changes with energy 
are present. Pulsed emission is detected up to ~ 24 keV. Note the enhancement 
within the main valley between the two peaks near phase 0.45, also visible in 
Fig. 7c of Gavriil & Kaspi (2002). 



monies are required to adequately fit the profiles in the various 
energy slices, mainly driven by the existence of the narrow, sig- 
nificant feature in the main valley between the peaks. This was 
also shown by Gavriil & Kaspi (2002). While the power of the 
second harmonic (2 maxima per cycle) dominates over that of 
the first up to ~ 10 keV (beyond which they become statistically 
equally strong) and the phase angle <£> 2 of the second harmonic 
remains more or less at the same position in phase, the first har- 
monic (broad component) shifts gradually to the right in phase 
with increasing energy (see Fig. 9). This energy dependent shift 
reflects the collapse at higher energies of the dominant pulse at 
low energies near phase 0.7. 

6.2. IE 2259+586 spectral characteristics 

In this section the time-averaged pulsed spectrum of IE 2259 + 
is derived based on the RXTE PCA & HEXTE observations 
given in Table 1 . The high-energy picture has been completed 
by comparing our results with the total emission spectrum (2-10 
keV) as determined by Patel et al. (2001) and with IBIS ISGRI 
(20-300 keV) and COMPTEL (0.75-30 MeV) upper limits (den 
Hartog et al. 2006). 

6.2. 1 . RXTE PCA/HEXTE pulsed spectrum 



E [keV] 

FIG. 9. — Phase angles as a function of energy for the first 3 harmonics 
used in the truncated Fourier series fit of the RXTE PCA pulse profiles of IE 
2259 + 586. Note that the position of the phase angle of the dominant second 
harmonic (labelled 2) is more or less constant (energy independent), while that 
of the first harmonic (labelled 1) decreases with increasing energy (shifts to 
the right). 



Employing the same counts-extraction technique as used above 
for the other AXPs, the pulsed excess counts from the PCA pro- 
files are converted into flux values taking into account the dif- 
ferent PCU exposures and sensitivities, and adopting an absorb- 
ing Hydrogen column of N H = (9.3 ±0.3) x 10 21 cm -2 (Patel et 
al. 2001). 

Assuming a simple absorbed power-law model the fit resulted 
in a poor x 2 of 26. 24 for 14-2 degrees of freedom (x 2 = 2.187), 
indicating an inappropriate fitting function (there is a ~ 1% 
probability that such a high value of x 2 is obtained at ran- 
dom assuming that the fit function is a good representation of 
the unkown parent function). An absorbed double power-law 
model, however, yielded a reasonable x 2 of 16.70 for 14-4 
degrees of freedom (~ 10% random probability). The improve- 
ment Ax 2 of 9.55 adding two additional fit parameters trans- 
lates to a ~ 3cr improvement adopting a maximum likelihood 
ratio test. Thus, the double power-law model provided a signif- 
icant improvement in describing the same spectral data. Adopt- 
ing this double-power-law model as underlying photon model 
spectrum, the pulsed flux values are shown in Fig. 10 as aqua 
colored filled circles, and the best fit double-power-law model 
as a dashed aqua colored line. 

The soft and hard power-law indices are T\ = 4.26 ± 0.01 and 
T2 = -1.02+0 [j, respectively, and the power-law-model com- 
ponents become equally strong at E- mt = 15.8 ±2.3 keV. Thus, 
also for this AXP we witness the onset of a dramatic hardening 
of the pulsed spectrum beyond ~ 10 keV, although confirma- 
^ggtion for energies beyond 20 keV is required. For HEXTE the 
source is too weak to be detected; 2cr upper limits (dark blue) 
are shown in Fig. 10. 

6.2.2. Discussion of spectra and pulsed fractions of IE 
2259 + 586 

The ISGRI flux upper limits to the total emission from IE 
2259 + 586 for energies above 20 keV (den Hartog et al. 2006) 
can be directly compared with the Chandra ACIS spectrum for 
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FIG. 10. — A vF v spectral representation of the total (2 keV - 30 MeV) 
and pulsed (2-30 keV) high-energy emission from IE 2259 + 586. The dashed 
aqua coloured line represents the double power-law fit to the PCA 2-30 keV 
pulsed spectrum. This fit suggests the onset of a hard spectral tail in the pulsed 
spectrum near ~ 15 keV. 

the 2-10 keV range. The black-body plus power-law model fit 
to this spectrum (Patel et al. 2001) is shown in Fig. 10 (dark- 
orange solid line). The double-power-law spectral fit suggests 
that again a minimum in luminosity is reached around 20 keV, 
where the extrapolation of the Chandra model fit crosses the 
pulsed PCA spectrum. Given the hard pulsed component sug- 
gested in these PCA data up to ~ 24 keV, we estimate that the 
source will show up at higher energies in an IBIS ISGRI mo- 
saic totalling > 4 Ms in effective on-axis exposure, what will 
be reached when adding in follow-up work already performed 
and/or scheduled observations. To complete the high-energy 
picture, the COMPTEL 0.75-30 MeV upper limits to the total 
emission (den Hartog et al. 2006) are also included in Fig. 10. 

Comparing in the overlapping energy range below 10 keV the 
PCA time-averaged pulsed fluxes with the Chandra total emis- 
sion model, and assuming no time variability, indicates a rather 
high value of > 43% for the pulsed fraction (pulsed/total emis- 
sion), in agreement with earlier estimates (Patel et al. 2001). 

7. IE 1048.1 -5937 

IE 1048. 1 - 5937 was discovered with Einstein at an angular 
separation of ~ 40' from r\ Car (Seward & Chlebowski 1982). 
During an observation on July 13, 1979 a sinusoidally shaped 
signal with a high pulsed fraction of 68 ±7% was detected at 
a period of ~ 6.44 s (Seward & Charles 1984; Seward et al. 
1986). The pulsations at a period of 6.4407(9) s were confirmed 
in EXOSAT ME observations (Smale et al. 1985; Seward et 
al. 1986), and the positional accuracy of IE 1048.1 -5937 was 
further tightened up down to a radius of about 10". 

Ginga observations combined with earlier measurements re- 
vealed a steady increase in spin period at a mean rate of P = 
(4.64 ± 1.1) x lO^syr" 1 (Corbet & Day 1990), implying a ro- 
tational energy loss much smaller than the lower limit on the 
X-ray luminosity. 

X-ray observations in the early nineties by ROS AT and ASCA 



(Mereghetti 1995; Corbet & Mihara 1997) indicated that the 
spin-down rate almost doubled after 1988, likely anticorrelated 
with the X-ray flux. 

BeppoSAX LECS / MECS data from long exposures in May 
1997 indicated the necessity of a two component spectral model, 
a power-law plus black body, to describe properly the measured 
X-ray (0.5-10 keV) spectrum (Oosterbroek et al. 1998). This 
was confirmed by Paul et al. (2000) using ASCA data. In the 
meantime, detailed timing studies using RXTE monitoring ob- 
servations performed in 1997-2000 (Kaspi et al. 2001) showed 
that significant deviations from simple spin-down exist, making 
phase-coherent timing impossible over time stretches longer 
than a few months. Inspite of these rotational irregularities, nei- 
ther pulse profile changes nor large pulsed flux variations were 
found. IE 1048. 1-5937 exhibited three X-ray bursts, on 2001 
Oct. 29, 2001 Nov. 14, and 2004 June 29, all caught during 
RXTE monitoring observations (Gavriil et al. 2002b, 2006). 

Precise X-ray imaging of IE 1048.1-5937 with Chandra 
(Wang & Chakrabarty 2002; Israel et al. 2002) provided a source 
position with sub-arcsecond accuracy. Within the 0."7 error 
circle only one faint near-IR source was found which clearly 
showed variability. A detailed variability study in the infrared 
and optical was undertaken by Durant & van Kerkwijk (2005). 
They established the variable nature in the infrared and optical 
and found a possible anticorrelation with the X-ray pulsed flux 
(see Gavriil & Kaspi 2004, for an extensive study of the X-ray 
pulsed flux variability and spin-down rate variations). At radio 
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FIG. 11.— RXTE PCA pulse profiles of IE 1048.1 -5937 for energies in the 
range 2.2-11.9 keV combining data collected between July 27, 1996 and Feb. 
24, 2004 (see Table 1). Two cycles are shown for clarity. 
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frequencies (21 cm) an expanding hydrogen shell centered on 
IE 1048.1-5937 was detected (Gaensler et al. 2005), which 
can be interpreted as a wind bubble blown by a 30-40 Mq star, 
likely the massive progenitor of IE 1048.1 -5937. 

Finally, the large effective area of XMM Newton, combined 
with accurate imaging, made it possible to perform in-depth 
(phase-resolved) spectral analyses not polluted by effects caused 
by the presence of strong nearby X-ray sources. XMM ob- 
served IE 1048. 1 - 5937 on three occasions and results are pre- 
sented by Tiengo et al. (2002); Mereghetti et al. (2004) and 
Tiengo et al. (2005). Comparing the three XMM observations 
revealed long-term flux and pulsed fraction variations in anti- 
correlation. The featureless spectral shape, however, remained 
more or less the same, and could be described by a combination 
of a power-law with photon index T ~ 2.7-3.5 and a blackbody 
with temperature kT ~ 0.63 keV. Phase-resolved spectroscopy 
clearly indicated that the spectrum is softer at pulse minimum 
and harder at pulse maximum with respect to the phase aver- 
aged spectrum. 

In the next section we present for the first time the hard X- 
ray/soft 7-ray timing and spectral characteristics of IE 1048. 1 - 
5937 above ~ 8 keV using a) all available public RXTE PCA 
and HEXTE data (see Table 1), b) INTEGRAL IBIS ISGRI 
data (see Table 2) and c) CGRO COMPTEL 0.75-30 MeV data. 

7.1. IE 1048.1 - 5937 timing characteristics 

In the timing analysis of the full set of PCA observations of 
IE 1048.1 -5937 listed in Table 1 according to the guidelines 
given in Sect. 3.1, we produced time-averaged pulse-phase dis- 
tributions for energies between ~ 2 and 35 keV. In the fold- 
ing/correlation process (see Sect. 3. 1) we used the ephemerides 
given by Kaspi et al. (2001). Significant pulsed emission has 
been detected up to ~ 12 keV. The PCA pulse profiles of IE 
1048.1-5937 are shown in Fig. 11 for the three following 
differential energy bands, 2.2-4.0, 4.0-7.9 and 7.9-11.9 keV. 
The detection of pulsed emission, with a non-uniformity sig- 
nificance of 6.8(7, between ~ 7.9 and 11.9 keV is reported for 
the first time. The single peaked pulse profiles, with a slightly 
steeper rise than fall, do not show morphology changes as a 
function of energy. HEXTE data do not show pulsed emission 
at energies above ~ 15 keV. 

7.2. IE 1048.1 - 5937 spectral characteristics 

7.2.1. RXTE PCA/HEXTE pulsed spectrum 

The derived pulsed excess counts have been converted to flux 
values assuming an absorbing column density of Nh = 1 .0 x 
10 22 cm" 2 (see Tiengo et al. 2005). It turned out that a single 
power-law model for the PCA band did not give an acceptable 
fit. Therefore we used for this AXP a black body plus power- 
law input spectrum. The optimized model (x 2 = 0.55 for 16-4 
d.o.f.; kT = 0.717(4) keV and T = 2.93(7)) is shown in Fig. 12 
together with the flux measurements, both aqua colored. It is in- 
teresting to note, that, while the photon index is in the expected 
range (Tiengo et al. 2005), the (time-averaged) black body tem- 
perature of the pulsed component is slightly higher than the val- 
ues obtained from the other X-ray instruments. This is very 
likely caused by the lack of sensitivity of the PCA for energies 
less than ~ 2.5 keV, thus poorly constraining the black body 
model at energies below its maximum. For HEXTE we could 
only derive upper limits for the pulsed emission, also shown in 
Fig. 12 as blue symbols. 
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FIG. 12. — A vF v representation of the spectrum of IE 1048.1-5937. 
Above ~ 12 keV the source is too weak to be detected so far. The IBIS ISGRI 
(purple) and CGRO COMPTEL (red) 2cr upper-limits are for the total emis- 
sion. Below 10 keV we also added spectral information for the pulsed (dotted) 
and total (solid) emission from two XMM Newton observations performed at 
different epochs (Tiengo et al. 2005), illustrating the variable nature both in 
normalization and shape of both components. 



7.2.2. INTEGRAL IBIS ISGRI and CGRO COMPTEL total 

spectrum 

IBIS ISGRI data (Core program - GPS or public) from ob- 
servations executed between revolution 30 and 217, satisfying 
our off-source pointing constraint of 14.°5 have been used to ob- 
tain spectral information on the total emission from IE 1048. 1 - 
5937 in the 20 - 300 keV energy range. The effective on-axis 
exposure time is, however, small ~ 250.8 ks, resulting in rather 
high upper-limits given the non-detections in any of the chosen 
broad energy bands. These 2a upper-limits have been included 
in Fig. 12 as purple symbols. Also, the analysis of full mis- 
sion COMPTEL data (4.8 Ms exposure in total) did not result 
in significant detections of IE 1048.1-5937. The 2er COMP- 
TEL upper-limits are shown as red symbols in Fig. 12. It is clear 
from the full high-energy spectral picture presented in Fig. 12 
that current spectral information above ^10 keV does not ex- 
clude the presence of a hard spectral component. Additional 
IBIS ISGRI data e.g. from a deep 2 Ms observation of the Ca- 
rina region, performed between INTEGRAL revolutions 192 
and 203, with IE 1048.1-5937 always in the FCFOV, will be 
very usefull to constrain the spectral properties of IE 1048.1 - 
5937 further in a future study. 

8. IE 1841-045 

IE 1841-045 is the first AXP for which surprisingly non- 
thermal pulsed X-ray/soft 7-ray emission was discovered (Kuiper 
et al. 2004); see this paper also for a summary of earlier obser- 
vational results on IE 1841-045). The pulsed spectrum above 
10 keV could be described, fitting RXTE PCA and HEXTE 
data, by a power-law model with photon index 0.94 ± 0.16 up 
to ~ 150 keV. The total spectrum of IE 1841-045 presented 
in Kuiper et al. (2004) was derived from HEXTE data using 
source ON/OFF rocking-mode observations. We realized later 
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Pulse Phase 

FIG. 13.— RXTE PCA pulse profiles of IE 1841-045 for the 2.1-10.3 and 
10.3-20.1 keV energy ranges (panels a and b, respectively), compared with the 
IBIS ISGRI 20-100 keV pulse profile (panel c). The dotted lines at phases 0.7 
and 1.05 serve as a guide to the eye for alignment comparison purposes. The 
20-100 keV IBIS ISGRI profile mimics the 10.3-20.1 keV RXTE PCA profile, 
suggesting a smooth prolongation above 20 keV of the 10.3-20.1 keV profile 
shape. 



that the HEXTE flux values were contaminated by a contribu- 
tion from the Galactic diffuse X-ray emission (see discussion in 
section 3.2). Better spectral information can be obtained with 
an imaging instrument like INTEGRAL IBIS ISGRI. Prelimi- 
nary total flux measurements with ISGRI were so far only re- 
ported by Molkov et al. (2004). For the present work signif- 
icantly more ISGRI data are available. Therefore, we present 
in this section a) a new spectrum of the total emission and b) 
performed now also a timing analysis for IE 1841-045 using 
ISGRI data. Furthermore, c) we repeated the timing/spectral 
RXTE PCA analysis using additional data compared to Kuiper 
et al. (2004) as well as improved event selection criteria (see 
Sect. 3.1). Finally, d) also for this source we analyzed COMP- 
TEL sky maps. 

8.1. IE 1841-045 pulse profiles for RXTE PCA and 
INTEGRAL IBIS ISGRI 

In our RXTE reanalysis for IE 1841-045, data from observa- 
tion 80098 (spread over 28-3-2003 - 17-2-2004; screened PCU- 
2 exposure ~ 69. 1 ks) were added to the set used in Kuiper et 
al. (2004) (now, total screened PCU-2 exposure: ~ 340.2 ks). 
We used in the folding/correlation process (see Sect. 3.1) the 
ephemerides given by Gotthelf et al. (2002) and derived in this 
work (see Table 3). The results are pulse profiles with improved 
statistics, confirming the variation of pulse shape with energy 
as presented by (Kuiper et al. 2004). As an example, Fig. 13 
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FIG. 14. — A vF v spectral representation of the high-energy emission from 
IE 1841-045 and SNR Kes 73. The aqua, blue and magenta data points/curves 
represent the (time averaged) pulsed emission component of the spectrum (2- 
300 keV) based on RXTE PCA/HEXTE and INTEGRAL IBIS ISGRI ob- 
servations. The black dashed line shows the best power-law model fit to the 
PCA, HEXTE and ISGRI pulsed flux values for energies above ~ 15 keV. The 
purple data points/curves reflect the IBIS ISGRI total flux measurements and 
fits (power-law, dashed; curved power-law, dashed-dotted and extrapolated to 
750 keV). The red 2<t upper limits from CGRO COMPTEL require a spectral 
break, somewhere between 140 and 750 keV. For completeness, below 10 keV 
also the total flux from IE 1841-045 plus SNR Kes 73 (yellow symbols; XMM 
Newton) and the total flux from IE 1841-045 (dark orange solid line; Chandra 
ACIS) are shown (see Kuiper et al. 2004, for more details). 



(panels a and b) shows RXTE PCA profiles for the 2.1-10.3 and 
10.3-20.1 keV energy bands. 

IBIS ISGRI data from observations with source off-axis an- 
gles less than < 14.°5 and taken during INTEGRAL revolutions 
49 - 123 have been processed to yield pulse phase distributions 
for differential energy bands between 20 and 300 keV. The to- 
tal screened effective on-axis exposure was ~ 808.8 ks. We 
derived a very accurate (RMS residual is 0.015 in phase; see 
Table 3 for the details) phase coherent ephemeris using RXTE 
PCA monitoring observations performed during March - De- 
cember 2003. This has been applied in the folding process of 
the selected barycentered ISGRI events. For the 20-100 keV 
energy range we obtained a Z\ of 3.6er. This ISGRI pulse pro- 
file is also shown in Fig. 13 (panel c). Its shape mimics that of 
the PCA 10.3-20.1 keV profile suggesting a stable pulse shape 
in the hard X-ray window. The statistics are still low, however, 
but additional ISGRI data from already performed and future 
scheduled observations will allow more detailed studies of the 
pulse shape in this energy range. 

8.2. IE 1841-045 spectral characteristics 

8.2.1. RXTE and INTEGRAL pulsed emission 

Starting with the new RXTE PCA phase distributions, intro- 
duced above, the updated IE 1841-045 pulsed spectrum, now 
with significant flux values up to ~ 28 keV, is shown in Fig. 14 
using aqua colored symbols. Also given are the earlier-reported 
pulsed-flux values from HEXTE (dark blue data points; see 
Kuiper et al. 2004), which clearly revealed the hard pulsed com- 
ponent up to ~ 150 keV. 



16 



Kuiper et al. 



Now we can also derive information on the pulsed spectrum 
from IBIS ISGRI. Phase distributions have been generated, fol- 
lowed by extraction of pulsed counts, for three logarithmically 
spaced energy bands (20-50, 50-120 and 120-300 keV), yield- 
ing only for the 50-120 keV band a significant pulsed-flux mea- 
surement. This ISGRI pulsed-flux value and two upper limits, 
shown in Fig. 14 as magenta symbols, are fully consistent with 
the hard spectrum measured by HEXTE. We fitted the new PCA 
pulsed spectrum for energies above ~ 15 keV together with the 
HEXTE and INTEGRAL pulsed flux data points with a sim- 
ple power-law model, yielding as best fit a photon index T of 
0.72 ±0.15, slightly harder than derived by Kuiper et al. (2004). 
This best fit model is indicated as a dashed black line in Fig. 14. 

8.2.2. INTEGRAL total emission 

In this work we derived the total emission spectrum of IE 
1841-045 in the 20-300 keV energy range from IBIS ISGRI 
mosaic maps generated for 7 differential energy bands over 
the 20 - 300 keV band, thus significantly improving upon the 
spectral information given by Molkov et al. (2004). We used a 
combination of public and Core program observations, all with 
source off -pointing angles less than 14.°5, and for which the de- 
tails 6 are listed in Table 2. The total screened effective on-axis 
exposure on IE 1841-045 in this combination was ~ 1,111 Ms. 
Significant emission has been detected up to ~ 140 keV (95- 
140 keV; 5.7 ± 0.9 mCrab). The derived flux values are shown 
in Fig. 14 as purple colored symbols. Over this 20-300 keV 
band we can not discriminate between a power-law or power- 
law model with energy dependent index to be the best model de- 
scribing the measurements. However, including in Fig. 14 the 
CGRO COMPTEL flux upper limits (red symbols) derived in 
this work, we notice that the best fitting power-law model with 
photon index V = 1 .32 ±0.11 does not extend into the MeV do- 
main. The spectrum must break somewhere between 140 and 
750 keV. The curved power-law model is consistent with the 
COMPTEL 2cr upper limits. Both models are shown in Fig. 14 
as purple dashed (power-law) and purple dashed-dotted (curved 
power-law) lines with the latter model extrapolated up to 750 
keV. 

8.2.3. Discussion of spectra and pulsed fractions of IE 
1841-045 

The new and better total spectrum of IE 1841-045, now de- 
rived using ISGRI data, is indeed lower, ~ 40 % between 20 
and 100 keV, than published earlier using HEXTE data (Kuiper 
et al. 2004), and is also slightly harder. It connects better to 
the total spectra measured at different epochs by XMM New- 
ton and Chandra ACIS below 10 keV, shown in Fig. 14 and 
adopted from Kuiper et al. (2004, and citations therein). 

For the pulsed spectrum, the power-law fit to the new PCA 
spectrum above ~ 15 keV combined with the new ISGRI mea- 
surements and the earlier published HEXTE flux values, con- 
firms the drastic up turn/hardening of the pulsed spectrum. Also, 
the extrapolation of this very hard spectrum has to remain un- 
der our COMPTEL upper limits, requiring an even more drastic 
bend/break compared to what is required for the total spectrum. 

The compilation of Chandra, XMM-Newton, RXTE and IN- 
TEGRAL spectra in Fig. 14, taken at very different epochs 
over many years, suggests that the hard X-ray emission of IE 

6 Data from Revs. 131,133-134 were ignored because these have been taken 
in a staring mode configuration giving rise to substantial systematical effects in 
the image deconvolution procedures. 



1841-045 is stable, i.e. this magnetar is most of the time in the 
same state of activity for energies above ^10 keV. The pulsed 
fraction is confirmed to be ~ 25% at 20 keV and <~ 100% for 
energies beyond ^100 keV. 

9. SUMMARY 

Exploiting the availability of archival data from RXTE mon- 
itoring observations and new, deep INTEGRAL exposures of 
AXPs, we were able to show that AXPs exhibit exceptionally 
hard spectra for energies above 10 keV. Of the sample of five 
AXPs studied in this work, three (the brightest at energies be- 
low lOkeV: IE 1841-045, 1RXS J1708-4009, 4U 0142 + 61) 
are shown to emit up to energies of <~ 150 keV. Of the two 
weaker sources, for one (IE 2259 + 586) an upturn of pulsed 
emission is suggested up to ~ 25 keV, and the other (IE 1048.1 - 
5937) is too weak to detect a possible similar upturn in the hard- 
X-ray range with presently available exposures. We regard this 
sufficient evidence to conclude that a persistent non-thermal 
component at energies above 10 keV is a common property of 
AXPs. 

In Table 4 we summarize for the five AXPs the spectral pa- 
rameters presented in this work, together with estimates for 
their distances and pulsar characteristics (period, period deriva- 
tive, deduced surface magnetic field strength at the pole and 
spin-down power). Most interesting and new are the results on 
the total and pulsed spectra, and luminosities above 10 keV. As 
noted above, the exceptionally hard pulsed spectra seem to ex- 
tend up to at least 150 keV, therefore we have given for all four 
AXPs for which the spectral up turn around 10 keV has been 
measured, the luminosity of the pulsed component integrating 
the best-fit power-law model over the full decade in energy from 
10 to 100 keV. For the luminosity of the pulsed component be- 
tween 1 and 10 keV, we list the power in the power-law com- 
ponent fitted to the pulsed spectra before the break, except for 
IE 1048.1 -5937, where the power has been derived from a 
black body plus power-law model. Thus for 4 AXPs we ig- 
nore the presence of a black-body component which is known 
to be present at lower energies as well. However, for 1RXS 
J1708-4009, 4U 0142 + 61, IE 2259 + 586 and IE 1841-045 
we could fit the pulsed PCA spectra down to about 2 keV well 
with a single power-law model. 

The total spectra above 20 keV appear also to be hard, but 
somewhat softer than the pulsed spectra. 

For IE 1841-045 and 4U 0142 + 61 the total emission be- 
comes consistent with 100% pulsed around ^100 keV, starting 
from a pulsed fraction at 10 keV of ~ 25% and ~ 10% for 
IE 1841-045 and 4U 0142 + 61, respectively. The exception is 
1RXS J 1708 -4009, for which the hard X-ray emission above 
10 keV is consistent with being 100% pulsed. Note in this con- 
text, however, the strong intensity variability reported for this 
AXP for energies below 10 keV, which might be coupled to 
variations in pulsed fraction. 

Table 4 immediately shows that the luminosities of the hard 
X-ray components (pulsed and total) exceed the available total 
spin-down powers by a few orders of magnitude, a conclusion 
drawn earlier for the total luminosities (black body plus power- 
law components) in the 1-10 keV band. 

Our compilations of total and pulsed flux values for energies 
above 10 keV for the three strongest AXPs, using data from 
different instruments and collected at different epochs, suggest 
that the high-energy component is stable, contrary to reports for 
the emissions below 10 keV. However, addional observations 
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Pulsar characteristics and high-energy spectral properties (pulsed/total) of the AXPs studied in this work. 
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are required (and will become available with INTEGRAL) to 
constrain this further. 

The time-averaged light curves for the different AXPs sug- 
gest that the pulse shapes for energies above 10 keV exhibit less 
structure, and vary less with energy than seen below 10 keV. 
However, the statistics are still insufficient to draw stringent 
conclusions. Phase-resolved spectroscopy on time-averaged pro- 
files will be performed when more INTEGRAL data can be 
added. 

10. DISCUSSION 

The discussion of AXPs in the context of magnetically pow- 
ered rotating neutron stars, was till recently focussed on under- 
standing the observational characteristics of the soft-spectrum 
component for energies below 10 keV. Indeed, evidence for 
hard X-ray or soft gamma-ray emission from AXPs was lack- 
ing, in contrast to the very luminous outbursts at soft gamma- 
ray energies characterizing SGRs. Our discovery of pulsed 
emission from AXPs up to at least ^150 keV adds a completely 
new non-thermal component requiring a steady mechanism for 
accelerating particles in magnetospheres of magnetars. Also 
this new quiescent emission component is far too luminous to 
be powered by rotational energy loss, as is evident in Table 4. 
Recently, INTEGRAL also found for the first time quiescent 
hard X-ray emission from a SGR (SGR 1806-20; Molkov et al. 
2005; Mereghetti et al. 2005). As for AXPs, its spectrum above 
10 keV could be fitted with a hard power-law model, photon 
index 1.6 ± 0.1, an index comparable to, but somewhat softer 
than that measured for the total AXP spectra in that energy win- 
dow, but significantly softer than the hard spectra reported in 
this work for the quiescent pulsed spectra of AXPs. In addi- 
tion, the quiescent SGR spectrum does not exhibit a black body 
component at energies below 10 keV. This underlines again, 
that AXPs and SGRs are (very) different manifestations of the 
magnetar scenario. 

In Table 4, we listed the high-energy characteristics for five 
AXPs studied in this work, out of a total of six persistent AXPs 
known to date. Although this sample is still very small, we 
verified whether there are already any apparent correlations be- 
tween the measured parameters. There is a hint for a correla- 
tion between the magnetic field strength B s and the luminosity 
^no-ioo °f tne P u l se d emission of the hard tail (10-100 keV) for 
four AXPs, as well as the total emission of the hard tail (now 
for three AXPs). However, IE 1048.1-5937 has a B s three 
times higher than 4U 0142 + 61, and is located at a similar dis- 
tance, but no hard X-ray emission has been detected, and its 
luminosity at ~ 10 keV is about three times lower than that of 
4U 0142 + 61. It is interesting to note that the spin down lumi- 
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FIG. 15. — A vFy spectral representation of the pulsed emission for the 4 
AXPs showing hard spectral tails. For comparison purposes, also the pulsed 
high-energy (~ 1 keV - 10 GeV) spectra of the young Crab (PSR B0531+21) 
and PSR B 1509-58 pulsars, and the middle-aged Vela pulsar (PSR B0833- 
045) are shown. Furthermore, the 3tr GLAST sensitivity (green dashed-dotted 
line) assuming an E 2 source for an all-sky survey duration of one year (> 100 
MeV) is shown. 



nosity of IE 1048.1 - 5937 is the highest of this sample, mainly 
due to its large P, and reaches almost 50% of the luminosity of 
the pulsed emission between 1 and 10 keV, a factor 10 - 100 
higher than for the other four AXPs. 

Surprisingly, the high-energy AXP spectra are very similar 
to those of "middle-aged" radio pulsars in the ~ 10-300 keV 
energy band. This can be seen in Fig. 15 in which the AXP 
spectra are shown together with the spectra of two young ra- 
dio pulsars, PSR B0531+21 (Crab) and PSR B 1509-58, and the 
middle-aged Vela pulsar (PSR B0833-45). These are the only 
radio pulsars which have been detected from soft X-rays up to 
the gamma-ray domain. The other middle-aged high-energy 
gamma-ray pulsars detected by CGRO EGRET (see e.g. the 
review by Thompson et al. 1997) are detected at high-energy 
gamma-rays (>100 MeV), and soft X-rays, but not in the in- 
termediate hard X-ray range, in which they are weaker than the 
Vela pulsar. The overall high-energy spectral shape of these 
middle-aged pulsars is expected to be similar to that of Vela. 
Like the AXPs, the Vela pulsar has a (pulsed) thermal black- 
body component (only visible below 2 keV, and not shown in 
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Fig. 15), and an extrapolation of its power-law shape spectrum 
above 10 keV to the IR, optical and radio domains is approx- 
imately in agreement with the measured flux values at these 
longer wavelengths. This is rather similar to the case of e.g. 
AXP 4U 0142 + 61, with the difference that none of the AXPs 
have so far been seen in the radio, and the X-ray and IR emis- 
sions exhibit a time variable behaviour, not seen for Vela. This 
similarity of the multiwavelengths spectra could suggest some 
similarity in production processes in the two different scenarios 
of rotation-powered pulsars and magnetars starting from very 
different sources of power. 

A first attempt, before our discovery of hard X-ray emission 
from AXPs, to model possible production of high-energy (X- 
ray and gamma-ray) emission in the magnetospheres of AXPs, 
applying a scenario proposed for high-energy production in the 
magnetospheres of radio pulsars, was made by Cheng & Zhang 
(2001). They modeled the production of high-energy (above 
100 MeV) gamma-radiation in outer magnetospheric gaps of 
AXPs. They argued that due to the strong field of a magne- 
tar, the gamma-ray emission rooted at the polar caps will be 
quenched. However, far away from the pulsar surface, i.e. in 
outer vacuum gaps, gamma radiation could be emitted because 
the local field will drop below the critical quantum limit. They 
predict that the NASA Gamma-ray Large Area Space Telescope 
GLAST, due for launch in 2007, might be able to detect this 
very-high-energy emission from AXPs (see Fig. 15 for the 3a 
GLAST sensitivity assuming an E~ 2 source spectrum for an all- 
sky survey duration of one year). However, their model calcu- 
lations do not reproduce the hard spectra and high X-ray lumi- 
nosities now discovered above 10 keV. 

The publication of the discovery of a hard spectral tail of 
IE 1841-045 by Kuiper et al. (2004) stimulated Thompson & 
Beloborodov (2005) to reconsider high-energy emission from 
magnetars due to an ultra strong magnetic field. A gradual re- 
lease of energy in the stellar magnetosphere is expected if it is 
twisted and a strong electric current is induced on the closed 
field lines. They considered two mechanisms of gamma-ray 
emission: 

(1) A thin surface layer of the star is heated by the down- 
ward beam of current carrying charges, which excite Lang- 
muir turbulance in the layer. Thus, a temperature kT of ~ 
100 keV can be reached, emitting bremsstrahlung photons up 
to this characteristic temperature. Interestingly, three of the 
AXPs, 1RXS J1708-4009, 4U 0142 + 61 and IE 1841-045, 
were detected up to these energies of <~ 150 keV, the shapes 
of their total and/or pulsed spectra being consistent with this 
mechanism with peak luminosities around the predicted energy. 
However, the 100% pulsed fraction around 20 keV for 1RXS 
J 1708 -4009 is difficult to reconcile with the surface emission, 
given the effects of general relativistic light bending. Future 
observations have to reveal whether the cut off in the spectra is 
also consistent with this interpretation. 

(2) Soft 7-rays are produced at a distance of ~ 100 km from 
the star surface in the magnetosphere, where the electron cy- 
clotron energy is in the keV range. A large electric field de- 
velops in this region in response to the outward drag force felt 
by the current-carrying electrons from the flux of keV photons 
leaving the star. Thompson & Beloborodov show that a seed 
photon injected in this region undergoes a runaway accelera- 
tion and upscatters keV photons above the threshold for pair 
creation. The resulting synchrotron spectrum can reach a peak 
at ~ 1 MeV. Our hard-X-ray spectra together with the upper 



limits from COMPTEL in the MeV domain, allow spectra with 
maximum luminosities close to the MeV range. But we should 
note that the flux values reported around ~ 100 keV are already 
at about the same level as the 2-er upper limits for the 0.75- 
3 MeV interval, suggesting a different shape, thus challenging 
this interpretation. 

An alternative quantum-electrodynamics model for the bursts 
and the quiescent non-thermal emission from AXPs as well as 
SGRs was proposed by Heyl & Hernquist (2005a,b). Their 
model is based on fast-mode breakdown, in which wave energy 
is dissipated into electron-positron pairs when the scale of these 
discontinuities becomes comparable to an electron Compton 
wavelength. They showed that under appropriate conditions, 
an electron-positron fireball would result, producing primarily 
the X-ray and soft-gamma-ray bursts. They also investigated 
(Heyl & Hernquist 2005b) the properties of non-thermal emis- 
sion in the absence of a fire ball. This quiescent, non-thermal 
emission associated with fast-mode breakdown may account 
for the observed non-thermal emission presented in this work 
for AXPs, as well as for the quiescent emission reported for 
SGR 1806-20 (Molkov et al. 2005; Mereghetti et al. 2005). In- 
deed, they succeeded in fitting ISGRI AXP spectra as well as 
unabsorbed optical data for 4U 0142 + 61 from Hulleman et al. 
(2000), and the SGR 1806-20 non-thermal spectrum above 10 
keV. Interestingly, they predict that the emission should extend 
beyond the observed INTEGRAL spectra without a break be- 
low 1 MeV. However, the combination of ISGRI spectra and 
COMPTEL upper limits presented in this work, seem to con- 
tradict this claim. They further state, that if the magnetars have 
significant optical excesses, such as 4U 0142 + 61, then the qui- 
escent emission from most of the AXPs discussed in this work 
and of SGR 1806-20 should be detectable with GLAST (see 
again Fig. 15 for the GLAST sensitivity). Further INTEGRAL 
observations, and ultimately GLAST observations can support 
or reject these claims. 

Over the last few years, two radio pulsars with periods and 
period time derivatives similar to those found for AXPs, thus 
also with high-magnetic-field strengths, are found: PSR J 1847- 
0130, which has the highest by far inferred surface dipolar mag- 
netic field (B = 9.4 x 10 13 G) among all known radio pulsars 
(McLaughlin et al. 2003) and PSR J1718-3718 (Hobbs et al. 
2004; Kaspi & McLaughlin 2005). The X-ray luminosities (or 
upper-limits to these) in the 2-10 keV range derived for these 
two radio pulsars are much lower than those of the standard 
AXP group. Above 10 keV, detection of hard X-rays could 
connect these two different types of neutron stars. PSR J 1847- 
0130 and PSR J1718-3718 both happened to be located near 
the centers of our deep IBIS ISGRI exposures on IE 1841-045 
and 1RXS J1708-4009, respectively. Therefore, we searched 
above 20 keV for hard X-ray signatures from these pulsars, but 
nothing was detected in any of the energy bands we have inves- 
tigated. The 2cr upper limit on the flux in the 20-30 keV band 
for PSR J1847-0130 is 0.4 mCrab 7 and that for PSR J1718- 
3718 is also 0.4 mCrab, but now in the 20-35 keV band. 

Therefore, why solitary rotating neutron stars with AXP-like 
timing parameters in some cases behave as "dull" radio pulsars 
and in other cases as "enigmatic" magnetars is still unclear. 
Different suggestions have been made (see e.g. McLaughlin et 
al. 2003), one of the possibilities being that high-field pulsars 
and AXPs have similar dipole magnetic fields, but AXPs have 
also quadrupole (or higher) components. Kaspi & McLaughlin 

7 IE 1841-045 is detected in the same image with a flux of ~ 1 .5 ± 0.2 mCrab 
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(2005) suggest a.o. the interesting possibility that the high-field 
radio pulsars, may one day emit transient magnetar-like emis- 
sion, and conversely that the transient AXPs might be more 
likely to exhibit radio pulsations. More theoretical work is re- 
quired, but certainly also more detailed observational results at 
all wavelengths to verify the model predictions. Ongoing obser- 
vations of AXPs with the wide-field-of-view INTEGRAL IBIS 
imager will allow us to contribute with deeper studies in the 
hard X-ray/soft 7-ray range. 

This research has made use of data obtained from the High 
Energy Astrophysics Science Archive Research Center (HEAS ARC) 
provided by NASA's Goddard Space Flight Center, and of data 
obtained through the INTEGRAL Science Data Centre (ISDC), 
Versoix, Switzerland. INTEGRAL is an ESA project with in- 
struments and science data centre funded by ESA member states 
(especially the PI countries: Denmark, France, Germany, Italy, 
Switzerland, Spain), Czech Republic and Poland, and with the 
participation of Russia and the USA. We have extensively used 
NASA's Astrophysics Data System (ADS). We thank Anton 
Klumper and Jaap Schuurmans for the maintenance of the OS A 
software at SRON and the INTEGRAL data import from the 
ISDC site. Nanda Rea is acknowledged for useful discussions 
on AXPs in general. Finally, we appreciate the willingness of 
Jacco Vink and Maurizio Falanga to make their INTEGRAL 
open time data of the Cassiopeia region (Revs. 261-269) di- 
rectly available for us. 
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